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Abstract. Quantum control is concerned with active manipulation of physical 
and chemical processes on the atomic and molecular scale. This work presents a 
perspective of progress in the field of control over quantum phenomena, tracing the 
evolution of theoretical concepts and experimental methods from early developments 
to the most recent advances. Among numerous theoretical insights and technological 
improvements that produced the present state-of-the-art in quantum control, there 
have been several breakthroughs of foremost importance. On the technology side, 
the current experimental successes would be impossible without the development of 
intense femtosecond laser sources and pulse shapers. On the theory side, the two 
most critical insights were (1) realizing that ultrafast atomic and molecular dynamics 
can be controlled via manipulation of quantum interferences and (2) understanding 
that optimally shaped ultrafast laser pulses are the most effective means for producing 
the desired quantum interference patterns in the controlled system. Finally, these 
theoretical and experimental advances were brought together by the crucial concept of 
adaptive feedback control, which is a laboratory procedure employing measurement- 
driven, closed-loop optimization to identify the best shapes of femtosecond laser control 
pulses for steering quantum dynamics towards the desired objective. Optimization 
in adaptive feedback control experiments is guided by a learning algorithm, with 
stochastic methods proving to be especially effective. Adaptive feedback control 
of quantum phenomena has found numerous applications in many areas of the 
physical and chemical sciences, and this paper reviews the extensive experiments. 
Other subjects discussed include quantum optimal control theory, quantum control 
landscapes, the role of theoretical control designs in experimental realizations, and 
real-time quantum feedback control. The paper concludes with a prospective of open 
research directions that are likely to attract significant attention in the future. 
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1. Introduction 

For many decades, physicists and chemists have employed various spectroscopic methods 
to carefully observe quantum systems on the atomic and molecular scale. The fascinating 
feature of quantum control is the ability to not just observe but actively manipulate 
the course of physical and chemical processes, thereby providing hitherto unattainable 
means to explore quantum dynamics. This remarkable capability along with a multitude 
of possible practical applications have attracted enormous attention to the field of control 
over quantum phenomena. This area of research has experienced extensive development 
during the last two decades and continues to grow rapidly. A notable feature of this 
development is the fruitful interplay between theoretical and experimental advances. 

Various theoretical and experimental aspects of quantum control have been 
reviewed in a number of articles and books [1-15]. This paper starts with a short 
review of historical developments as a basis for evaluating the current status of the field 
and forecasting future directions of research. We try to identify important trends, follow 
their evolution from the past through the present, and cautiously project them into the 
future. This paper is not intended to be a complete review of quantum control, but 
rather a perspective and prospective on the field. 

In section 2, we discuss the historical evolution of relevant key ideas from the first 
attempts to use monochromatic laser fields for selective excitation of molecular bonds, 
through the inception of the crucial concept of control via manipulation of quantum 
interferences, and to the emergence of advanced contemporary methods that employ 
specially tailored ultrafast laser pulses to control quantum dynamics of a wide variety 
of physical and chemical systems in a precise and effective manner. After this historical 
summary, we review in more detail the recent progress in the field, focusing on significant 
theoretical concepts, experimental methods, and practical advances that have shaped the 
development of quantum control during the last decade. Section 3 is devoted to quantum 
optimal control theory (QOCT), which is currently the leading theoretical approach for 
identifying the structure of controls (e.g., the shape of laser pulses) that enable attaining 
the quantum dynamical objective in the best possible way. We present the formalism 
of QOCT (i.e., the types of objective functionals used in various problems and methods 
employed to search for optimal controls), consider the issues of controllability and 
existence of optimal control solutions, survey applications, and discuss the advantages 
and limitations of this approach. In section 4, we review the theory of quantum 
control landscapes, which provides a basis to analyze the complexity of finding optimal 
solutions. Topics discussed in that section include the landscape topology (i.e., the 
characterization of critical points), optimality conditions for control solutions, Pareto 
optimality for multi-objective control, homotopy trajectory control methods, and the 
practical implications of control landscape analysis. The important theoretical advances 
in the field of quantum control have laid the foundation for the fascinating discoveries 
occurring in laboratories where closed-loop optimizations guided by learning algorithms 
alter quantum dynamics of real physical and chemical systems in dramatic and often 
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unexpected way. Section 5, which constitutes a very significant portion of this paper, is 
devoted to laboratory implementations of adaptive feedback control (AFC) of quantum 
phenomena. We review numerous AFC experiments that have been performed during 
the last decade in areas ranging from photochemistry to quantum information sciences. 
These experimental studies (most of which employ shaped femtosecond laser pulses) 
clearly demonstrate the capability of AFC to manipulate dynamics of a broad variety 
of quantum systems and explore the underlying physical mechanisms. The role of 
theoretical control designs in experimental realizations is discussed in section 6. In 
particular, we emphasize the importance of theoretical studies for the feasibility analysis 
of quantum control experiments. Section 7 presents concepts and potential applications 
of real-time feedback control (RTFC). Both measurement-based and coherent types of 
RTFC are described, along with current technological obstacles limiting more extensive 
use of these approaches in the laboratory. Future directions of quantum control are 
considered in section 8, including important unsolved problems and some emerging new 
trends and applications. Finally, concluding remarks are given in section 9. 

2. Early developments of quantum control 

The historical origins of quantum control lie in early attempts to use lasers for 
manipulation of chemical reactions, in particular, selective breaking of bonds in 
molecules. Lasers, with their tight frequency control and high intensity, were considered 
ideal for the role of molecular-scale 'scissors' to precisely cut an identified bond, without 
damage to others. In the 1960s, when the remarkable characteristics of lasers were 
initially realized, it was thought that transforming this dream into reality would be 
relatively simple. These hopes were based on intuitive, appealing logic. The procedure 
involved tuning the monochromatic laser radiation to the characteristic frequency 
of a particular chemical bond in a molecule. It was suggested that the energy of 
the laser would naturally be absorbed in a selective way, causing excitation and, 
ultimately, breakage of the targeted bond. Numerous attempts were made in the 1970s 
to implement this idea [46-48]. However, it was soon realized that intramolecular 
vibrational redistribution of the deposited energy rapidly dissipates the initial local 
excitation and thus generally prevents selective bond breaking [49-51]. This process 
effectively increases the rovibrational temperature in the molecule in the same manner 
as incoherent heating does, often resulting in breakage of the weakest bond(s), which is 
usually not the target of interest. 

2.1. Control via two-pathway quantum interference 

Several important steps towards modern quantum control were made in the late 1980s. 
Brumer and Shapiro [52-55] identified the role of quantum interference in optical control 
of molecular systems. They proposed to use two monochromatic laser beams with 
commensurate frequencies and tunable intensities and phases for creating quantum 
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interference between two reaction pathways. The theoretical analysis showed that 
by tuning the phase difference between the two laser fields it would be possible to 
control branching ratios of molecular reactions [56-58]. The method of two-pathway 
quantum interference can be also used for controlling population transfer between bound 
states [59,60] (in this case, the number of photons absorbed along two pathways often 
must be either all even or all odd to ensure that the wave functions excited by the two 
lasers have the same parity; most commonly, one- and three-photon excitations were 
considered) . 

The principle of coherent control via two-pathway quantum interference was 
demonstrated during the 1990s in a number of experiments, including control of 
population transfer in bound-to-bound transitions in atoms and molecules [59- 
64], control of energy and angular distributions of photoionized electrons [65-68] 
and photodissociation products [69] in bound-to-continuum transitions, control of 
cross-sections of photochemical reactions [70-72], and control of photocurrents in 
semiconductors [73, 74]. However, practical applications of this method are limited 
by a number of factors. In particular, it is quite difficult in practice to match 
excitation rates along the two pathways, either because one of the absorption cross- 
sections is very small or because other competing processes intervene. Another practical 
limitation, characteristic of experiments in optically dense media, is undesirable phase 
and amplitude locking of the two laser fields [75]. Due to these factors and other technical 
issues (e.g., imperfect focusing and alignment of the two laser beams), modulation depths 
achieved in two-pathway interference experiments were modest: typically, about 25- 
50% for control of population transfer between bound states [60-62, 64] (the highest 
reported value was about 75% in one experiment [63]), and about 15-25% for control 
of dissociation and ionization branching ratios in molecules [70, 71]. Two-pathway 
interference control is a nascent form of full multi-pathway control offered by operating 
with broad-bandwidth optimally shaped pulses. 

2.2. Pump-dump control 

In the 1980s, Tannor, Kosloff, and Rice [76, 77] proposed a method for selectively 
controlling intramolecular reactions by using two successive femtosecond laser pulses 
with a tunable time delay between them. The first laser pulse (the "pump" ) generates 
a vibrational wave packet on an electronically excited potential-energy surface of the 
molecule. After the initial excitation, the wave packet evolves freely until the second 
laser pulse (the "dump" ) transfers some of the population back to the ground potential- 
energy surface into the desired reaction channel. Reaction selectivity is achieved by 
using the time delay between the two laser pulses to control the location at which 
the excited wave packet is dumped to the ground potential-energy surface [7, 11]. For 
example, it may be possible to use this method to move the ground-state wave-function 
beyond a barrier obstructing the target reaction channel. In some cases, the second 
pulse transfers the population to an electronic state other than the ground state (e.g.. 
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to a higher excited state) in a pump-repump scheme. 

The feasibihty of the pump-dump control method was demonstrated in a number 
of experiments [78-82]. The pump-dump scheme can be also used as a time- resolved 
spectroscopy technique to explore transient molecular states and thus obtain new 
information about the dynamics of the molecule at various stages of a reaction [83-90] . 
In pump-dump control experiments, the system dynamics often can be explained in 
the time domain in a simple and intuitive way to provide a satisfactory qualitative 
interpretation of the control mechanism. The pump-dump method gained considerable 
popularity [7,11,30] due to its capabilities to control and investigate molecular dynamics. 
However, the employment of transform-limited laser pulses significantly restricts the 
effectiveness of this technique as a practical control tool. More effective control of the 
wave-packet dynamics and, consequently, higher reaction selectivity can be achieved 
by optimally shaping one or both of the pulses. For example, even a chirp of the 
pump pulse may improve the effectiveness of control by producing more localized wave 
packets (the use of pulse chirping will be discussed in section 2.6 in more detail). Recent 
experimental applications of the pump-dump scheme with shaped laser pulses (optimized 
using adaptive methods) will be discussed in section 5. 

2.3. Control via stimulated Raman adiabatic passage 

In the late 1980s, Bergmann and collaborators [91-94] demonstrated a very efficient 
adiabatic method for population transfer between discrete quantum states in atoms or 
molecules. In this approach known as stimulated Raman adiabatic passage (STIRAP), 
two time-delayed laser pulses (typically, of nanosecond duration) are applied to a three- 
level A-type configuration to achieve complete population transfer between the two lower 
levels via the intermediate upper level. Interestingly, the pulse sequence employed in 
the STIRAP method is counter- intuitive, i.e., the Stokes laser pulse that couples the 
intermediate and final states precedes (but overlaps) the pump laser pulse that couples 
the initial and intermediate states. The laser electric fields should be sufficiently strong 
to generate many cycles of Rabi oscillations. The laser-induced coherence between the 
quantum states is controlled by tuning the time delay, so that the transient population 
in the intermediate state remains almost zero, thus avoiding losses by radiative decay. 
Detailed reviews of STIRAP and related adiabatic passage techniques can be found 
in [10, 16]. While the efficiency of the STIRAP method, under appropriate conditions, 
is very high, its applicability is restricted to control of population transfer between a 
few discrete states as arise in atoms and small (diatomic and triatomic) molecules. In 
larger polyatomic molecules, the very high density of levels generally prevents successful 
adiabatic passage [10,16]. 

2.4. Control via wave-packet interferometry 

Another two-pulse approach for control of population transfer between bound states 
employs Ramsey interference of optically excited wave packets [95,96]. In this method, 
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referred to as wave-packet interferometry (WPI) [12], two time-delayed laser pulses 
excite an atomic, molecular, or quantum-dot transition, resulting in two wave packets 
on an excited state. Quantum interference between the two coherent wave packets 
can be controlled by tuning the time delay between the laser pulses. For control of 
population transfer, constructive or destructive interference between the excited wave 
packets gives rise to larger or smaller excited-state population, respectively. The same 
control mechanism is also applicable to other problems such as control of atomic 
radial wave-functions and control of molecular alignment. WPI was demonstrated 
with Rydberg [97-99] and fine-structure [100, 101] wave packets in atoms, vibrational 
[102-106] and rotational [107] wave packets in molecules, and exciton fine-structure wave 
packets in semiconductor quantum dots [108, 109] (for a detailed review of coherent 
control applications of WPI, see [12]; the use of WPI for molecular state reconstruction 
is reviewed in [110]). Once again, much more effective manipulation of quantum 
interferences is possible in this control scheme when shaped laser pulses are used instead 
of transform-limited ones (see section 5 for details). 

2.5. Introduction of Q OCT 

Although the control approaches discussed in sections 2.1-2.4 were initially perceived 
as quite different, it is now clear that on a fundamental level all of them employ the 
mechanism of quantum interference induced by control laser fields. A common feature 
of these methods is that they generally attempt to manipulate the evolution of quantum 
systems by controlling just one parameter: the phase difference between two laser fields 
in control via two-pathway quantum interference; the time delay between two laser 
pulses in pump-dump control, STIRAP, and WPI. While single-parameter control may 
be relatively effective in some simple systems, more complex systems and applications 
require more flexible and capable control resources. The single-parameter control 
schemes have been unifled and generalized by the concept of control with specially 
tailored ultrashort laser pulses. Rabitz and co-workers [111-113] and others [114,115] 
suggested that it would be possible to steer the quantum evolution to a desired product 
channel by speciflcally designing and tailoring the time- dependent electric fleld of the 
laser pulse to the characteristics of the system. Specifically, QOCT may be used to 
design laser pulse shapes which are best suited for achieving the desired goal [111-122]. 
An optimally shaped laser pulse typically has a complex form, both temporally and 
spectrally. The phases and amplitudes of different frequency components are optimized 
to excite an interference pattern amongst distinct quantum pathways, to best achieve the 
desired dynamics. The first optimal fields for quantum control were computed by Shi, 
Woody, and Rabitz [111] who showed that the amplitudes of the interfering vibrational 
modes of a laser-driven molecule could add up constructively in a given bond. We will 
review QOCT and its applications in more detail in section 3 (for earlier reviews of 
QOCT, see [7,11,14,38,40]). 
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2.6. Control with linearly chirped pulses 

Laser pulse-shaping technology rapidly developed during the early 1990s [4, 5, 12]. 
However, the capabilities of pulse shaping were not fully exploited in quantum control 
until the first experimental demonstrations of adaptive feedback control (AFC) in 
1997-1998 [123, 124]. Initially, ultrashort laser pulses with time- varying photon 
frequencies were used to tune just the linear chirp, which represents an increase 
or decrease of the instantaneous frequency as a function of time under the pulse 
envelope. § Linearly chirped femtosecond laser pulses were successfully applied for 
control of various atomic and molecular processes, including control of vibrational 
wave packets [125-131], control of population transfer between atomic states [132-134] 
and between molecular vibrational levels [135-137] via "ladder-climbing" processes, 
control of electronic excitations in molecules [138-142], selective excitation of vibrational 
modes in coherent anti-Stokes Raman scattering (CARS) [143], improvement of the 
resolution of CARS spectroscopy [144,145], and control of photoelectron spectra [146] 
and transitions through multiple highly excited states [117] in strong-field ionization 
of atoms. In particular, when the emission and absorption bands of a molecule 
strongly overlap, pulses with negative and positive chirp excite vibrational modes 
predominately in the ground and excited electronic states, respectively [125,129-131]. 
Chirped pulses can be also used to control the localization of vibrational wave packets in 
diatomic molecules, with the negative and positive chirp increasing and decreasing the 
localization, respectively [126-128]. Based on this effect, pump pulses with negative 
chirp were used to enhance selectivity in pump-dump control of photodissociation 
reactions [128]. Recently, the localization effect of negatively chirped pulses was used 
to protect vibrational wave packets against rotationally-induced decoherence [148]. Due 
to their effectiveness in various applications, chirped laser pulses are widely used in 
quantum control. However, by the end of the 1990s, many experimenters realized that 
more sophisticated pulse shapes, beyond just linear chirp, provide a much more powerful 
and flexible tool for control of quantum phenomena in complex physical and chemical 
systems. Femtosecond pulse-shaping technology is utilized to the fullest extent in AFC 
experiments where laser pulses are optimally tailored to meet the needs of complex 
quantum dynamics objectives [13,15,18,19,22,25-29,31,37,41]. The enormous growth 
of this field during the last decade is reviewed in section 5. 

2. 1. Control via non-resonant dynamic Stark effect 

Optimal control of quantum phenomena in atoms and molecules usually operates at laser 
intensities sufficient to be in the non-perturbative regime. Thus, controlled dynamics 
will naturally utilize the dynamic Stark shift amongst other available physical processes 
in order to reach the target. In a recent quantum control development, Stolow and co- 
workers proposed and experimentally demonstrated manipulation of molecular processes 

§ The instantaneous frequency uj(t) of a linearly chirped pulse with a carrier frequency wo is given at 
time t by uj{t) = ujq + 2bt, where b is the chirp parameter that can be negative or positive. 
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exclusively employing the non- resonant dynamic Stark effect (NRDSE) [149-152]. In 
this approach, a quantum system is controlled by an infrared laser pulse in the 
intermediate field-strength regime (non-perturbative but non-ionizing). Laser frequency 
and intensity are chosen to eliminate the complex competing processes (e.g., multiphoton 
resonances and strong-field ionization), so that only the NRDSE contributes to the 
control mechanism. By utilizing Raman coupling, control via NRDSE reversibly 
modifies the effective Hamiltonian during system evolution, thus making it possible 
to affect the course of intramolecular dynamic processes. For example, a suitably 
timed infrared laser pulse can act as a "photonic catalyst" by reversibly modifying 
potential energy barriers during a chemical reaction without inducing any real electronic 
transitions [152]. Control via NRDSE was successfully applied to create field- free 
"switched" wave packets (which can be employed, e.g., for molecular axis alignment) 
[149, 150] and modify branching ratios in non-adiabatic molecular photodissociation 
[151,152]. 

2. 8. Control of nuclear spins with radiofrequency fields 

One of the earliest examples of coherent control of quantum dynamics is manipulation 
of nuclear spin ensembles using radiofrequency (RF) fields [153]. The main application 
of nuclear magnetic resonance (NMR) control techniques is high-resolution spectroscopy 
of polyatomic molecules (e.g., protein structure determination) [154-157]. While control 
of an isolated spin by a time-dependent magnetic field is a simple quantum problem, 
in reality, NMR spectroscopy of molecules containing tens or even hundreds of nuclei 
involves many complex issues such as the effect of interactions between the spins, thermal 
relaxation, instrumental noise, and influence of the solvent. Therefore, modern NMR 
spectroscopy often employs thousands of precisely sequenced and phase-modulated 
pulses. Among important NMR control techniques are composite pulses, refocusing, 
and pulse shaping. In particular, the use of shaped RF pulses in NMR makes it possible 
to improve the frequency selectivity, suppress the solvent contribution, simplify high- 
resolution spectra, and reduce the size and duration of experiments [158]. In recent years, 
NMR became an important testbed for developing control methods for applications in 
quantum information sciences [159-163]. In order to perform fault-tolerant quantum 
computations, the system dynamics must be controlled with an unprecedented level of 
precision, which requires even more sophisticated designs of control pulses than in high- 
resolution spectroscopy. In particular, QOCT was recently applied to identify optimal 
sequences of RF pulses for operation of NMR quantum information processors [163, 164]. 

3. Quantum optimal control theory 

The most comprehensive means for coherently controlling the evolution of a quantum 
system (e.g., a molecule) undergoing a complex dynamical process is through the 
coordinated interaction between the system and the electromagnetic field whose spectral 
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content and temporal profile may be continuously altered throughout the process. For a 
specified control objective, and with restrictions imposed by many possible constraints, 
the time-dependent field required to manipulate the system in a desired way can be 
designed using QOCT [14,38,40]. This general formulation encompasses both the weak 
and strong field limits and incorporates as special cases the one-parameter methods such 
as control via two-pathway quantum interference and pump-dump control. 

3.1. Description of controlled quantum dynamics 

Optimal control theory has an extensive history in traditional engineering applications 
[165, 166], but control of quantum phenomena imposes special features. Consider the 
time evolution of a controlled quantum system, 

±Uit) = -^[Ho-eit).f^]Uit), UiO) = L (1) 

Here, U(t) is the unitary evolution operator of the system at time t, I is the identity 
operator, Hq is the free Hamiltonian, fi is the dipole operator, and e(t) is the control 
function at time t. For the control function, we will use the notation e{-) G K (where 
IK is the space of locally bounded, sufficiently smooth, square integrable functions of 
time defined on some interval [0,T], with T fixed). Equation (1) adequately describes 
the coherent quantum dynamics of a molecular system interacting with a laser electric 
field in the dipole approximation or a spin system interacting with a time-dependent 
magnetic field. We will consider finite-level quantum systems and denote the dimension 
of the system's Hilbert space as N.\\ 

For quantum systems undergoing environment ally- induced decoherence, there are 
many dynamical models depending on the character of the system-environment coupling. 
If the system and environment are initially uncoupled, the evolution of the system's 
reduced density matrix p from t = to t is described by a completely positive, trace 
preserving map: p{t) = $(po) where po = p(0). This map (which is often called the 
Kraus map) can be expressed using the operator-sum representation (OSR) [167]: 

n 

p{t) = ^{p,) = Y,K,{t)p,K]{t), (2) 
i=i 

where Kj are Kraus operators (A^ x complex matrices), which satisfy the condition 
Yl^j=i ^j^j ~ -^A^' denotes the identity operator on the Hilbert space of 

dimension A^. There exist infinitely many OSRs (with different sets of Kraus operators) 
for the same Kraus map. Various types of quantum master equations can be derived 
from the Kraus OSR under additional assumptions [168, 169]. In particular, for the 

II If one is not concerned with the physically irrelevant global phase of the evolution operator, the 
control problem can be restricted to the Hamiltonian represented by a traceless Hermitian matrix, and 
U £ SU(A^), where SU(A^) is the special unitary group. 
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generic class of Markovian environments, the dynamics of an open quantum system can 
be described using the quantum master equation of the Lindblad type [168-170]: 



^p(t) = -^[Ho - eit) ■ /i,p(t)] + 5^ 7. L,p{t)Ll - \L\Up{t) - \p{t)L\L 

i=l '- 



(3) 



where 7^ are non-negative constants and are Lindblad operators (A^ x complex 
matrices) that represent the non-unitary effect of coupling to the environment. For a 
closed system, (3) reduces to the von Neumann equation, p{t) = —{i/h)[Ho—e{t)-p, p{t)]. 
For simplicity, we will formulate QOCT below using the Schrodinger equation (1) 
for the unitary evolution operator; analogous formulations using the von Neumann 
equation or the Lindblad master equation for the density matrix are available in the 
literature [120,171-174]. 



3.2. Control objective functionals 

The general class of control objective functionals (cost functionals of the Bolza type) 
can be written as 

Jm-),ei-)] = F(f/(T)) + r GiU{t),e{t))dt, (4) 

Jo 

where F is a continuously differentiable function on U(A^), and G is a continuously 
differentiable function on U(A^) x M. The optimal control problem may be stated as the 
search for 

Jopt = max J[f/(-),e(-)], (5) 

e(-) 

subject to the dynamical constraint (1). If only the term G{U{t),e{t))dt is present, 
the cost functional is said to be of the Lagrange type, whereas if only the term F{U (T)) 
is present, the functional is said to be of the Mayer type [166]. Three classes of problems 
corresponding to different choices of F{U{T)) have received the most attention in the 
quantum control community to date: (i) evolution-operator control, (ii) state control, 
and (iii) observable control. 

For evolution-operator control, the goal is to generate U{T) such that it is as close 
as possible to the target unitary transformation W. The Mayer-type cost functional in 
this case can be generally expressed as 

F,iU{T)) = l-\\W-U{T)l (6) 

where || ■ || is an appropriate normalized matrix norm, i.e., Fi{U(T)) is maximized when 
the distance between U{T) and W is minimized. This type of objective is common in 
quantum computing applications [i"")], where Fi{U{T)) is the fidelity of a quantum 
gate [176, 177]. One frequently used form of the objective functional Fi{U(T)) is 
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obtained utilizing the squared Hilbert-Schmidt norm [178] in (6) with an appropriate 
normahzation (i.e., ||X|| = {2N)-^Tt{X^X)) [179-181]: 

F,{UiT)) = j^^TT[W^U{T)]. (7) 

Other forms of the objective functional, which employ different matrix norms in (6), 
are possible as well [181-184]. For example, a modification of (7), Fi{U{T)) = 

|Tr [iy^f/(T)] I , which is independent of the global phase of f/(T), can be used. Note 
that Fi{U(T)) is independent of the initial state, as the quantum gate must produce 
the same unitary transformation for any input state of the qubit system [175]. 

For state control, the goal is to transform the initial state po into a final state p(T) = 
U{T)pqU\T) that is as close as possible to the target state pf. The corresponding 
Mayer-type cost functional is 

F^{U{T)) = 1 - \\U{T)p,U\T) - pfll, (8) 

where || ■ || is an appropriate normalized matrix norm (e.g., the Hilbert-Schmidt norm 
can be used) [185-188]. 

For observable control, the goal is typically to maximize the expectation value of a 
target quantum observable G (represented by a Hermitian operator). The corresponding 
Mayer-type cost functional is [116,189-192] 

F^{U{T)) = Tr [U{T)p^U\T)Q] . (9) 

An important special case is state-transition control (also known as population transfer 
control), for which po = ^"^^ © = \ipi){ipi\, where \ipi) and iV'f) are eigenstates 

of the free Hamiltonian in this case, the objective functional (9) has the form 
F^{U{T)) = Pi^f = \{ip{\U{T)\ipi)\'^, which is the probability of transition (i.e., the 
population transfer yield) between the energy levels of the quantum system [193, 194]. 
In many chemical and physical applications of quantum control, absolute yields are 
not known, and therefore maximizing the expectation value of an observable (e.g., the 
population transfer yield) is a more appropriate laboratory goal than minimizing the 
distance to a target expectation value. 

Also, in quantum control experiments (see section 5), measuring the expectation 
value of an observable is much easier than estimating the quantum state or evolution 
operator. Existing methods of quantum-state and evolution-operator estimation rely 
on tomographic techniques [195-202] which are extremely expensive in terms of the 
number of required measurements (e.g., in quantum computing applications, standard 
methods of state and process tomography require numbers of measurements that scale 
exponentially in the number of qubits [175,202-206]). Therefore, virtually all quantum 
control experiments so far have used observable control with objective functionals of the 
form (9). For example, in an AFC experiment [148], in which the goal was to maximize 
the degree of coherence, the expectation value of an observable representing the degree 
of quantum state localization was used as a coherence "surrogate," instead of state 
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purity or von Neumann entropy which are nonhnear functions of the density matrix and 
hence would require state estimation. Nevertheless, future laboratory applications of 
quantum control, in particular in the field of quantum information sciences, will require 
evolution-operator control and state control, with the use of objective functionals of 
the types (6) and (8), respectively, together with novel state and process estimation 
methods [205-214]. 

Recently, attention has turned to problems requiring simultaneous maximization 
of several control objectives [215-219]. In the framework of QOCT, these optimization 
problems are sometimes handled through the use of a weighted-sum objective functional, 
such as [216, 217] 

n 

F^{U{T)) = ^ afcTr [f/(T)pof/^(T)0fc] , (10) 

k=l 

which extends (9) to multiple quantum observables. Also, general methods of multi- 
objective optimization [220-222] have been recently applied to various quantum control 
problems [217-219]. 

Another common goal in quantum control is to maximize a Lagrange-type cost 
functional subject to a constraint on U{T) [223,224]. For example, this type of control 
problem can be formulated as follows: 

max / G{e{t)) dt, subject to F{U{T)) = x, (H) 
Jo 

where F{U(T)) is the Mayer-type cost functional for evolution-operator, state, or 
observable control (as described above), and x is a constant that corresponds to the 
target value of F. Often, the goal is to minimize the total field fiuence, in which case 
G{e{t)) = -ie^(t) is used. 

3.3. Controllability of quantum systems 

One of the fundamental issues of quantum control is to assess the system's controllability. 
A quantum system is controllable in a set of configurations, S = {(}, if for any pair of 
configurations Ci)C2 £ <S there exists a time-dependent control e{-) that can drive the 
system from the initial configuration (i to the final configuration ^2 in a finite time T. 
Here, the notion of configuration means either the state of the system p, the expectation 
value of an observable Tr(/?9), the evolution operator U, or the Kraus map $, depending 
on the specific control problem. Controllability of closed quantum systems with unitary 
dynamics has been well studied [225-244]. Controllability analysis was also extended to 
open quantum systems [245-252]. 

Controllability is determined by the equation of motion as well as properties of 
the Hamiltonian. For a closed quantum system with unitary dynamics (1), evolution- 
operator controllability implies that for any unitary operator W there exists a finite 
time T and a control e{-), such that W = U{T), where U{T) is the solution of (1). 
For an A^-level closed system, a necessary and sufficient condition for evolution-operator 
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controllability is that the dynamical Lie group & of the system (i.e., the Lie group 
generated by the system's Hamiltonian) be U(A^) (or SU(A^) for a traceless Hamiltonian) 
[235-237]. 

Unitary evolution preserves the spectrum of the quantum state (i.e., the eigenvalues 
of the density matrix). All density matrices that have the same eigenvalues form a 
set of unitarily equivalent states (e.g., the set of all pure states). Therefore, under 
unitary evolution, a quantum system can be state controllable only within a set of 
unitarily equivalent states [235,236]. Density-matrix controllability means that for any 
pair of unitarily equivalent density matrices pi and p2 there exists a control e{-) that 
drives pi into p2 (in a finite time). It has been shown [235-237] that density-matrix 
controllability is equivalent to evolution-operator controllability. For specific classes of 
density matrices, the requirements for controllability are weaker [229-231]. For example, 
pure-state controllability requires that the system's dynamical Lie group is transitive 
on the sphere §^^~^. For infinite-level quantum systems evolving on non-compact Lie 
groups, such as those arising in quantum optics, the conditions for controllability are 
more stringent [243,244,251,253]. 

3.4- Searching for optimal controls 

To solve for optimal controls that maximize an objective functional (of the types 
discussed in section 3.2), it is convenient to define a functional J that explicitly 
incorporates the dynamical constraint (1): 

J[Ui-), 0(-), ^(O] = FiUiT)) + \ r GiUit),eit))dt 

Jo 

-23f?^ Tr|0^(t) 

Here, A is a scalar weight and an auxiliary operator 0(t) is a Lagrange multiplier 
employed to enforce satisfaction of Eq. (1). 

Various modifications of the objective functional (12) are possible, for example, 
QOCT can be formulated for open systems with non-unitary dynamics [120, 171-174, 
254-258]. Modified objective functionals can also comprise additional spectral and 
fluence constraints on the control field [259,260], take into account nonlinear interactions 
with the control field [261,262], deal with time-dependent and time-averaged targets 
[256,263-265], and include the final time as a free control parameter [266,267]. It is 
also possible to formulate QOCT with time minimization as a control goal (time optimal 
control) [268-271]. As we mentioned earlier, QOCT can be also extended to incorporate 
optimization of multiple objectives [215-219]. 

A necessary condition for a solution of the optimization problem (5) subject to the 
dynamical constraint (1) is that the first-order functional derivatives of J with respect to 
U{-), 0(-), and e{-) are equal to zero. Correspondingly, optimal controls can be obtained 
by solving the resulting Euler-Lagrange equations. Equivalently, optimal controls can be 
derived through application of the Pontryagin maximum principle (PMP) [165,223,224]. 



e{t) ■ p) U{t) 



dt. 



(12) 
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Satisfaction of the first-order conditions following from the PMP is a necessary but not 
sufficient condition for optimality of a control e{-). So-called Legendre conditions on the 
Hessian, which depend on the type of cost, are also required for optimality [165,166]. 

An important issue is the existence of optimal control fields (i.e., maxima of 
the objective functional) for realistic situations that involve practical constraints on 
the applied laser fields. It is important to distinguish between the existence of an 
optimal control field and controllability; in the former field is designed, subject 

to particular constraints, that guides the evolution of the system towards a specified 
target until a maximum of the objective functional is reached, while in the latter 
case, the exact coincidence between the attained evolution operator (or state) and 
the target evolution operator (or state) is sought. The existence of optimal controls 
for quantum systems was analyzed in a number of works. Peirce et al. [11-] proved 
the existence of optimal solutions for state control in a spatially bounded quantum 
system that necessarily has spatially localized states and a discrete spectrum. Zhao and 
Rice [272] extended this analysis to a system with both discrete and continuous states 
and proved the existence of optimal controls over the evolution in the subspace of discrete 
states. Demiralp and Rabitz [185] showed that, in general, there is a denumerable 
infinity of solutions to a particular class of well-posed quantum control problems; the 
solutions can be ordered in quality according to the achieved optimal value of the 
objective functional. The existence of multiple control solutions has important practical 
consequences, suggesting that there may be broad latitude in the laboratory, even under 
strict experimental restrictions, for finding successful controls for well-posed quantum 
objectives. The existence and properties of critical points (including global extrema) 
of objective functionals for various types of quantum control problems were further 
explored using the analysis of control landscapes [179-181,189-194,273] (see section 4). 

A number of optimization algorithms were adapted or specially developed for 
use in QOCT, including the conjugate gradient search method [111], the Krotov 
method [177,274,275], monotonically convergent algorithms [276-282], non-iterative 
algorithms [283], the gradient ascent pulse engineering (GRAPE) algorithm [164], a 
hybrid local/global algorithm [256], and homotopy-based methods [284-286]. Faster 
convergence of iterative QOCT algorithms was demonstrated using "mixing" strategies 
[287]. Also, the employment of propagation toolkits [288-290] greatly increases the 
efficiency of numerical optimizations and allows for fast combinatorial optimization [291]. 
Detailed discussions of the QOCT formalism and algorithms are available in the 
literature [11,14,38,40]. 

3.5. An example of QOCT applied to a molecular system 

In order to illustrate optimal control of molecules, we consider an instructive example. 
In one of the pioneering QOCT studies, Kosloff et al. [11 1] considered two electronic 
states (ground and excited) of a model molecular system, with the wave function (in 



Control of quantum phenomena: Past, present, and future 



15 



the coordinate representation) of the form^ 

*(M) = ww)=(*j:;;|), (13) 

where ipg and ipe are the projections of the wave function on the ground and excited state, 
respectively. The time evolution of the wave function is determined by the Schrodinger 
equation: 

d_ ( M^,t) \( i/e(r) ifgc(r,t) A ( ^e(r,t) A 
dt[^,{r,t)) \Hl{T,t) H,{t) )\Mr,t))' ^ ^ 

where Hiijc) = p^/(2m) + Vi(r) {i = g, e), p is the momentum operator, V^(r) and 
Ve{r) are the adiabatic potential energy surfaces for the ground and excited state, 
respectively. The off-diagonal term Hge{r,t) represents the field-induced coupling 
between the molecular states: 

i/ge(r,t) = -^ge(r)£(t), (15) 

where /Uge(r) is the electric dipole operator and e(t) is the time-dependent electric field 
of the control laser pulse applied to the molecule. 

The goal is to control a dissociation reaction in the presence of two distinct exit 
channels on the ground potential energy surface. The corresponding objective functional 
(including the dynamical constraint) is given by 

J = {i;iT)\PmT)) - A l^'it) dt - 2^j\x{t)\ + W)) dt, (16) 

The first term in (16) represents the main control goal, where P is the projection 
operator on the state corresponding to the target exit channel (i.e., the part of the 
wave function which is beyond the target saddle point on the ground-state surface and 
is characterized by the outgoing momentum); the second term is used to manage the 
fluence of the control field, with A being a scalar weight factor; the third term includes 
an auxiliary state |x(^)) that is a Lagrange multiplier employed to enforce satisfaction 
of the Schrodinger equation [H is the 2x2 Hamiltonian matrix defined by (14)). In 
order to find the control field that maximizes the objective, the first-order functional 
derivatives of J with respect to x{')^ ^{') set to zero, producing the following 

Euler-Lagrange equations: 

ih^^\m) = H\m), \m) = \i^o), (17) 

i^^lx(t)) = H\x{t)), |x(T)) = P|^(T)), (18) 

% For the sake of notation consistency, the control problem is presented here slightly differently than 
in the original work [lf4]. 
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An initial guess is selected for the control field (e.g., a pair of transform-limited pulses 
with a time delay between them, as in pump-dump control), and equations (17)- 
(19) are solved using an appropriate algorithm, as discussed in section 3.4 above. 
This optimization procedure identifies a shaped control field eopt(") that maximizes 
photoinduced molecular dissociation into the target channel. Successful application of 
QOCT to this model molecular system [111] demonstrated the benefits of optimally 
tailoring the time-dependent laser field to achieve the desired dynamic outcome. 

3.6. Applications of QOCT 

Originally, QOCT was developed to design optimal fields for manipulation of molecular 
systems [111-122] and has been applied to a myriad of problems (e.g., rotational, 
vibrational, electronic, reactive, and other processes) [11, 14, 40]. Some recent 
applications include, for example, control of molecular isomerization [292-295], control 
of electron ring currents in chiral aromatic molecules [296], control of current flow 
patterns through molecular wires [297], and control of heterogeneous electron transfer 
from surface attached molecules into semiconductor band states [298]. Beyond 
molecules, QOCT has been applied to various physical objectives including, for 
example, control of electron states in semiconductor quantum structures [299-301], 
control of atom transport in optical lattices [302], control of Bose-Einstein condensate 
transport in magnetic microtraps [303], control of a transition of ultracold atoms from 
the superfluid phase to a Mott insulator state [-JOi], control of coherent population 
transfer in superconducting quantum interference devices [305], and control of the 
local electromagnetic response of nanostructured materials [306]. Recent interest has 
rapidly grown in applications of QOCT to the field of quantum information sciences, 
including optimal protection of quantum systems against decoherence [188,257,307-317], 
optimal operation of quantum gates in closed systems [176,177,224,253,318-332] and 
in open systems (i.e., in the presence of decoherence) [184,333-348], optimal generation 
of entanglement [266,267,347,349-351], and optimal (i.e., maximum-rate) transfer of 
quantum information [352]. In a recent experiment with trapped ion qubits, shaped 
pulses designed using QOCT were applied to enact single-qubit gates with enhanced 
robustness to noise in the control field [353]. Optimal control methods were also applied 
to the problem of storage and retrieval of photonic states in atomic media, including 
both theoretical optimization [354-356] and experimental tests [357-359]. 

3. 7. Advantages and limitations of QOCT 

An advantage of QOCT relative to the laboratory method of AFC (to be discussed 
in detail in section 5) is that the former can be used to optimize a well defined 
objective functional of virtually any form, while the latter relies on information obtained 
from measurements and thus is best suited to optimize expectation values of directly 
measurable observables. In numerical optimizations, there is practically no difference 
in effort between computing the expectation value of an observable, the density matrix. 
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or the evolution operator. In the laboratory, however, it is much more difficult to 
estimate a quantum state or, even more so, the evolution operator, than to measure the 
expectation value of an observable. Moreover, state estimation error increases rapidly 
with the Hilbert-space dimension [198,360]. The very large number of measurements 
required for accurate quantum state/process tomography [202-205] renders (at least, 
presently) the use of adaptive laboratory methods for state/evolution-operator control 
rather impractical (although not impossible). 

QOCT is often used to explore new quantum phenomena in relatively simple models 
to gain physical insights. The realization of quantum control is ultimately performed 
in the laboratory. In this context QOCT fits into what is called open-loop control. 
Generally, in open- loop control, a theoretical control design (e.g., obtained by using 
QOCT or another theoretical method) is implemented in the laboratory with the actual 
system. Unfortunately, there are not many problems for which theoretical control 
designs are directly applicable in the laboratory. QOCT is most useful when detailed 
knowledge of the system's Hamiltonian is available. Moreover, for open quantum 
systems, it is essential to know the details of the system-environment interaction. 
Therefore, the practical applicability of QOCT in the context of open-loop control is 
limited to very simple systems, i.e., mostly to cases when a small number of degrees 
of freedom can be controlled separately from the remainder of the system. This may 
be possible when the controlled subsystem has characteristic frequencies well separated 
from those of other transitions, and/or evolves on a time scale which is very different 
from that of the rest of the larger system. A well known example of such a separately 
controllable subsystem is nuclear spins in a molecule, which can be very well controlled 
using RF fields without disturbing rotational, vibrational, and electronic degrees of 
freedom. Another example is a subset of several discrete levels in an atom or diatomic 
molecule, the transitions between which can be controlled in a very precise way without 
any significant leakage of population to other states. However, for a majority of 
interesting physical and chemical phenomena, controlled systems are too complex and/or 
too strongly coupled to other degrees of freedom. For such complex systems, the 
accuracy of control designs obtained using model-based QOCT is usually inadequate, 
and hence laboratory AFC is generally the preferred strategy. In these situations, QOCT 
may be more useful for feasibility analysis and exploration of control mechanisms, as 
basic features of the controlled dynamics can be identified in many cases even using 
relatively rough models. 

4. Quantum control landscapes 

An important practical goal of quantum control is the discovery of optimal solutions 
for manipulating quantum phenomena. Early studies [112,185,272] described conditions 
under which optimal solutions exist, but did not explore the complexity of finding them. 
Underlying the search for optimal controls is the landscape which specifies the physical 
objective as a function of the control variables. Analysis of quantum control landscapes 
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[39] can not only establish the existence of optimal control solutions and determine their 
types (e.g., global versus local maxima and true maxima versus saddle points), but also 
deal with establishing necessary conditions for convergence of optimization algorithms 
to global maxima along with bounds on the scaling of convergence effort. Surprisingly, 
these properties are independent of details of a particular Hamiltonian (provided that 
the system is controllable), which makes the results of landscape analysis applicable 
across a wide range of controlled quantum phenomena. 



4.1. Control landscape definition and critical points 

Properties of the search space associated with Mayer-type cost functionals play a 
fundamental role in the ability to identify optimal controls. To characterize these 
properties, it is convenient to express the cost functional in a form where the dynamical 
constraints are implicitly satisfied. Consider a control problem with a fixed target time 
T for a closed quantum system with unitary evolution. Denote by Vr '■ b{-) h- U{T) 
the endpoint map from the space of control functions to the space of unitary evolution 
operators, induced by the Schrodinger equation (1), so that U(T) = Vt(£^(-))- A Mayer- 
type cost functional F{U{T)) itself describes a map F from the space of evolution 
operators to the space of real-valued costs. Thus the composition of these maps, 
J = F o Vt '■ K — > M, is a map from the space of control functions to the space 
of real-valued costs. This map generates the functional J[£:(-)] = -^(^r(^(')))- We will 
refer to the functional ^[£^(-)] the control landscape. The optimal control problem 
may then be expressed as the unconstrained search for 

Jopt = maxJ[e(-)]. (20) 

The topology of the control landscape (i.e., the character of its critical points, including 
local and global extrema) determines whether local search algorithms will converge to 
globally optimal solutions to the control problem [361]. Studies of quantum control 
landscape topology are presently an active research area [39, 179-181, 189, 191-194,233, 
362]. 

The critical points (extrema) of the landscape are controls, at which the first-order 
functional derivative of J[£^(-)] with respect to the control field is zero for all time, i.e., 

^^0. w.io.n. (21) 

The critical manifold Ai of the control landscape is the collection of all critical points: 

M = {e{-) I 5J/5e{t) =0, Vt G [0,T]} . (22) 

A central concept in landscape topology is the classification of a critical point as regular 
or singular [363,364]. Most generally, a critical point of ^[£^(-)] is regular if the map 
Vt is locally surjective in its vicinity, i.e., if for any local increment 6U{T) in the 
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evolution operator there exists an increment fe(-) in the control function such that 
Vrie^-) + Se{-)) = VT{e{-)) + 6U(T). This condition is equivalent to requiring that 
the elements /ijj(t) of the time-dependent dipole-operator matrix (in the Heisenberg 
picture) form a set of A^^ linearly independent functions of time [181]. In its turn, this 
condition is satisfied for all non-constant admissible controls if and only if the quantum 
system is evolution-operator controllable [181,226]. Note that for landscapes of some 
particular physical objectives the conditions for regularity of the critical points can be 
less stringent. For example, in the important special case of state-transition control, a 
critical point is regular if the matrix elements fiij(t) contain a set of just 2N — 1 linearly 
independent functions of time. This condition is satisfied for all non-constant admissible 
controls if and only if the quantum system is pure-state controllable (which is a weaker 
condition than evolution-operator controllability, as discussed in section 3.3). 

A critical point of J[£(-)] is singular if the map Vr is not locally surjective in the 
point's vicinity. Using the chain rule, one obtains: 



where VF(C/(r)) is the gradient of F at f/(T), 5U{T)/5e{t) is the first-order functional 
derivative of U{T) with respect to the control field, and {A, B) = Tt{A'''B) is the Hilbert- 
Schmidt inner product. From (23), if a critical point of J is regular, VF(U{T)) must 
be zero. A critical point is called kinematic if VF{U{T)) = and non-kinematic if 
VF{U{T)) 7^ 0. Thus, all regular critical points are kinematic. A singular critical point 
may be either kinematic or non-kinematic; in the latter case, 6J/6e{t) = whereas 
VF{U(T)) 7^ [364]. On quantum control landscapes, the measure of regular critical 
points appears to be much greater than that of singular ones [364]. Therefore attention 
has been focused on the characterization of regular critical points, and several important 
results have been obtained [39]. Nevertheless, singular critical points on quantum 
control landscapes have been recently studied theoretically [364] and demonstrated 
experimentally [365]. 

The condition for kinematic critical points, VF{U(T)) = 0, can be cast in an 
explicit form for various types of quantum control problems. For evolution-operator 
control with the objective functional J = Fi(f/(T)) of (7), this condition becomes 



i.e., W^U(T) is required to be a Hermitian operator. It was shown [179,180] that 
this condition implies W'^UlT) = y^(— © I]y^m)Y, where Y is an arbitrary unitary 
transformation and m = 0,1, . . . , N. There are A^ + 1 distinct critical submanifolds 
labeled by m, with corresponding critical values of J given by = 1 — {2m /N). 
The global optima corresponding to m = and m = N (with Jq = 1 and 
Jtv = —1, respectively) are isolated points, while local extrema corresponding to 
m = 1,2, . . . , N — 1 are smooth, compact, Grassmannian submanifolds embedded in 
U(A^). It can be shown that all regular local extrema are saddle-point regions [180]. 




(23) 



[179, 180] 



W^UiT) = U^{T)W, 



(24) 
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For observable control with the objective functional J = F3{U{T)) of (9), the 
condition for a kinematic critical point becomes [189, 193,233] 

[f/(T)pof/^(T),e]=0, (25) 

i.e., the density matrix at the final time is required to commute with the target 
observable operator. This condition was studied in the context of optimization of 
Lagrange-type cost functionals with an endpoint constraint [362, 366, 367] as well as 
in the context of regular critical points for Mayer- type cost functionals [189,233]. Let 
R and S denote unitary matrices that diagonalize po and 6, respectively, and define 
U{T) = S'^U{T)R. The condition (25) that p{T) and 6 commute is equivalent to the 
condition that the matrix U{T) is in the double coset Ai^ of some permutation matrix 
P. [191]: 

f/(T) eM^ = U(n)P^U(m). (26) 

Here, U(n) is the product group U(r2i) x ■■■ x U(r2f,), where U(n;) corresponds to 
the Ith eigenvalue of po with n;-fold degeneracy, and U(m) is the product group 
U(mi) X ■■■ X U(ms), where U(m;) corresponds to the /th eigenvalue of B with m;- 
fold degeneracy. Thus, each critical submanifold Ain corresponds to a particular choice 
of the permutation it. All permutations on indices form the symmetric group &n, 
and the entire critical manifold is given by A/f = U7re6]v "^^^ structure of A4 
depends on any degeneracies in the spectra of po and 6. When both po and 9 are 
fully nondegenerate, then U(n) = U(m) = [U(l)]^, and Ai consists of A^! disjoint A^- 
dimensional tori, labeled by the permutation matrices. The occurrence of degeneracies 
in the spectra of po and O will merge two or more tori together, thereby reducing the 
number of disjoint critical submanifolds and increasing their dimensions [191]. 



4-2. Optimality of control solutions 

Satisfaction of the condition (21) for a critical point is a necessary but not sufficient 
condition for optimality of a control [166,223]. For Mayer-type cost functionals, a 
sufficient condition for optimality is negative semidefiniteness of the Hessian of J, which 
is defined as 

The characteristics of critical points (in particular, the presence or absence of local 
optima) are important for the convergence properties of search algorithms [39]. To 
classify critical points as global maxima and minima, local maxima and minima, and 
saddle points, one examines the second-order variation in J for an arbitrary control 
variation Se{-), which for Mayer- type functionals can be written as 

52 J = Qp(SUiT), 5U{T)) + (VF(f/(T)), 5^U{T)), (28) 

where 5U{T) and 6'^U{T) are the first- and second-order variations, respectively, of 
U{T) caused by a control variation Se{-), and Qp is the Hessian quadratic form of 
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F{U). Assuming that the critical point e{-) is regular, one obtains: 

= Qp{5U{T),5U{T)). (29) 

Explicit expressions for the Hessian and/or Hessian quadratic form were obtained for 
evolution-operator control [179-181] and observable control [190,361]. 

The optimality of regular critical points can be determined by inspecting the 
number of positive, negative and null eigenvalues of the Hessian (or, equivalently, the 
coefficients of the Hessian quadratic form when written in a diagonal basis). An issue 
of special interest is to determine whether any of the regular critical points are local 
maxima (frequently referred to as local traps due to their ability to halt searches guided 
by gradient algorithms before reaching the global maximum). Detailed analyses for 
evolution-operator control and observable control reveal [39,179-181,190,361] that all 
regular optima are global and the remainder of regular critical points (i.e., except for 
the global maximum and global minimum) are saddles. This discovery means that no 
local traps exist in the control landscapes of controllable closed quantum systems. The 
same result was also obtained for observable-control landscapes of controllable open 
quantum systems with Kraus-map dynamics [273]. Due attention still needs to be given 
to consideration of singular critical points, although numerical evidence suggests that 
their effect on optimization is likely insignificant [36 1]. 

4-3. Pareto optimality for multi-objective control 

Many practical quantum control problems seek to optimize multiple, often competing, 
objectives. In such situations the usual notion of optimality is replaced by that of 
Pareto optimality. The Pareto front of a multi-objective control problem is the set 
of all controls such that all other controls have a lower value for at least one of the 
objectives [220-222]. The analysis of the Pareto front reveals the nature of confiicts 
and tradeoffs between different control objectives. The structure of the landscape for 
multi-observable control is of interest and follows directly from that of single-observable 
control [217]. Of particular relevance to many chemical and physical applications is the 
problem of simultaneous maximization of the expectation values of multiple observables. 
Such simultaneous maximization is possible if the intersection f]^ Ai^j^^^^ (where Ai^/^^^^ 
is the maximum submanifold for the kth observable) is nonempty and a point in the 
intersection can be reached under some control e{-)] in this regard, the dimension of 
the intersection manifold flfc-^i™'^^^ been analyzed [218]. It has been shown that 
the common QOCT technique of running many independent maximizations of a cost 
functional like (10) (using different weight coefficients {a/c}) is incapable of sampling 
many regions of the Pareto front ['218]. Alternative methods for Pareto front sampling 
are discussed further below. 
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4.4- Landscape exploration via homotopy trajectory control 

The absence of local traps in landscapes for observable control and evolution- 
operator control with Mayer-type cost functionals has important implications for the 
design of optimization algorithms. Many practical applications require algorithms 
capable of searching quantum control landscapes for optimal solutions that satisfy 
additional criteria, such as minimization of the field fluence or maximization of 
the robustness to laser noise. So-called homotopy trajectory control algorithms (in 
particular, diffeomorphic modulation under observable-response-preserving homotopy, 
or D-MORPH) [284-286] can follow paths to the global maximum of a Mayer-type 
cost functional, exploiting the trap-free nature of the control landscape, while locally 
optimizing auxiliary costs. The essential prerequisite for successful use of these 
algorithms is the existence of a connected path between the initial and target controls. 
Homotopy trajectory control is closely related to the notion of a level set which is defined 
as the collection of controls that all produce the same value of the cost functional J. 
Theoretical analysis [39, 285, 286] predicts that for controllable quantum systems each 
level set is a continuous manifold. A homotopy trajectory algorithm is able to move on 
such a manifold exploring different control solutions that result in the same value of the 
cost functional, but may differ in other properties (e.g., the field fluence or robustness). 
A version of the D-MORPH algorithm was also developed for evolution-operator control 
of closed quantum systems; it was able to identify optimal controls generating a target 
unitary transformation up to machine precision [324]. 

Homotopy trajectory algorithms are also very useful for exploring quantum control 
landscapes for multiple objectives. For example, in order to track paths in the space of 
expectation values of multiple observables while locally minimizing a Lagrange- type cost, 
multi-observable trajectory control algorithms were developed [217]. Such algorithms 
are generally applicable to the treatment of multi-objective quantum control problems 
(Pareto quantum optimal control) [218]. They can traverse the Pareto front to identify 
admissible tradeoffs in optimization of multiple control objectives (e.g., maximization 
of multiple observable expectation values). This method can continuously sample the 
Pareto front during the course of one optimization run [218] and thus can be more 
efficient than the use of standard QOCT with cost functionals of the form (10). Also, 
the D-MORPH algorithm was recently extended to handle optimal control problems 
involving multiple quantum systems and multiple objectives [219]. 

4-5. Practical importance of control landscape analysis 

The absence of local traps in control landscapes of controllable quantum systems has 
very important implications for the feasibility of AFC experiments (see section 5). 
The relationship between the quantum control landscape structure and optimization 
complexity of algorithms used in AFC has been the subject of recent theoretical 
analyses [39,368-370]. Results of these studies support the vast empirical evidence [361] 
indicating that the favorable landscape topology strongly correlates with fast mean 
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convergence times to the global optimum. The trap-free control landscape topology also 
ensures convergence of gradient-based optimization algorithms to the global maximum. 
These algorithms can be used to search for optimal solutions to a variety of quantum 
control problems. In addition to theoretical studies (mostly using QOCT), gradient 
algorithms are also applicable in quantum control experiments [371], provided that 
measurement of the gradient is sufficiently robust to laser and detection noise. The use of 
deterministic algorithms in AFC experiments is discussed in more detail in section 5.19. 

4-6. Experimental observation of quantum control landscapes 

Significant efforts have been recently devoted to experimental observation of quantum 
control landscapes, aiming both at testing the predictions of the theoretical analysis and 
at obtaining a better understanding of control mechanisms. Roslund et al. [372] observed 
quantum control level sets for maximization of non-resonant two-photon absorption in a 
molecule and second harmonic generation (SHG) in a nonlinear crystal and found them 
to be continuous manifolds (closed surfaces) in the control landscape. A diverse family 
of control mechanisms was encountered, as each of the multiple control fields forming 
a level set preserves the observable value by exciting a distinct pattern of constructive 
and destructive quantum interferences. 

WoUenhaupt, Baumert, and co-workers [373,374] used parameterized pulse shapes 
to reduce the dimensionality of the optimization problem (maximization of the 
Autler-Townes contrast in strong-field ionization of potassium atoms) and observed 
the corresponding two-dimensional quantum control landscape. In order to better 
understand the performance of AFC, the evolution of different optimization procedures 
was visualized by means of trajectories on the surface of the measured control landscape. 
Marquetand et al. [375] observed a two-dimensional quantum control landscape (for 
maximization of the retinal photoisomerization yield in bacteriorhodopsin) and used it 
to elucidate the properties of molecular wave-packet evolution on an excited potential 
energy surface. 

The theoretical analysis of control landscape topology has been carried out with 
no constraints placed on the controls (see section 4.1). A main conclusion from 
these studies is the inherent lack of local traps on quantum control landscapes under 
normal circumstances. Recently, Roslund and Rabitz [376] experimentally demonstrated 
the trap-free monotonic character of control landscapes for optimization of frequency 
unfiltered and filtered SHG. For unfiltered SHG, the landscape was randomly sampled 
and interpolation of the data was found to be devoid of traps up to the level of data noise. 
In the case of narrow-band-filtered SHG, trajectories taken on the landscape revealed 
the absence of traps, although a rich local structure was observed on the landscape in 
this case. Despite the inherent trap-free nature of the landscapes, significant constraints 
on the controls can distort and/or isolate portions of the erstwhile trap- free landscape 
to produce apparent (i.e., false) traps [376]. Such artificial structure arising from the 
forced sampling of the landscape has been seen in some experimental studies [373-375], 
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in which the number of control variables was purposely reduced. 
5. Adaptive feedback control in the laboratory 

There are important differences between quantum control theory and its experimental 
implementation. Control solutions obtained in theoretical studies strongly depend on 
the employed model Hamiltonian. However, for real systems controlled in the laboratory, 
the Hamiltonians usually are not known well (except for the simplest cases), and 
the Hamiltonians for the system-environment coupling are known to an even lesser 
degree. An additional difficulty is the computational complexity of accurately solving 
the optimal control equations for realistic polyatomic molecules. Another important 
difference between control theory and experiment arises from the difficulty of reliably 
implementing theoretical control designs in the laboratory, due to instrumental noise 
and other limitations. As a result, optimal theoretical control designs generally will 
not be optimal in the laboratory. Notwithstanding these comments, control simulations 
continue to be very valuable, and they even set forth the logic leading to practical 
laboratory control as explained below. 

A crucial step towards selective laser control of physical and chemical phenomena 
on the quantum scale was the introduction of AFC (also referred to as closed-loop 
laboratory control or learning control). AFC was proposed and theoretically grounded 
by Judson and Rabitz in their paper "Teaching lasers to control molecules" in 1992 [377] . 
In AFC, a loop is closed in the laboratory, with results of measurements on the quantum 
system used to evaluate the success of the applied control and to refine it, until the 
control objective is reached as best as possible. At each cycle of the loop, the external 
control (e.g., a shaped laser pulse) is applied to the system (e.g., an ensemble of 
molecules). The signal (e.g., the yield of a particular reaction product or population 
in a target state) is detected and fed back to the learning algorithm (e.g., a genetic 
algorithm). The algorithm evaluates each control based on its measured outcome with 
respect to a predefined control goal, and searches through the space of available controls 
to move towards an optimal solution. 

While AFC can be simulated on the computer [377-390], the important advantage 
of this approach lies in its ability to be directly implemented in the laboratory. Most 
importantly, the optimization is performed in the laboratory with the actual system, and 
thus is independent of any model. As a result, the AFC method works remarkably well 
for systems even of high complexity, including, for example, large polyatomic molecules 
in the liquid phase, for which only very rough models are available. Second, there is 
no need to measure the laser field in AFC, because any systematic characterization 
of the control "knobs" (such as pulse shaper parameters) is sufficient. This set of 
control "knobs" determined by the experimental apparatus defines the parameter space 
searched by the learning algorithm for an optimal laser shape. This procedure naturally 
incorporates any laboratory constraints on the control laser fields. Third, optimal 
controls identified in AFC are characterized by a natural degree of robustness to 
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instrumental noise, since non-robust solutions will be rejected by the algorithm. Fourth, 
in AFC, it is possible to operate at a high-duty cycle of hundreds or even thousands 
of experiments per second, by exploiting (i) the conceptual advantage of the evolving 
quantum system solving its own Schrodinger equation in the fastest possible fashion and 
(ii) the technological advantage of high-repetition-rate pulsed laser systems under full 
automation. Fifth, in AFC, a new quantum ensemble (e.g., a new molecular sample) 
is used in each cycle of the loop, which completely avoids the issue of back action 
exerted by the measurement process on a quantum system. Thus AFC is technologically 
distinct from measurement-based RTFC [391-396], in which the same quantum system 
is manipulated until the final target objective is reached and for which measurement 
back action is an important effect that needs to be taken into account (see section 7). 

5.1. Femtosecond pulse-shaping technology 

The majority of current AFC experiments employ shaped ultrafast laser pulses. In 
such experiments, one usually starts with a random or nearly random selection of trial 
shaped pulses of length ~ 10~^^ s or less. The pulses are shaped by modulating the 
phases and/or amplitudes of the spatially resolved spectral components, for example, 
by means of a liquid crystal modulator (LCM), acousto-optic modulator (AOM), or 
a micromechanical mirror array (MMA). The experiments employ fully automated 
computer control of the pulse shapes guided by a learning algorithm. The shaped laser 
pulses produced by this method can be viewed as "photonic reagents," which interact 
with matter at the atomic or molecular scale to facilitate desired controlled outcomes 
of various physical and chemical phenomena. 

Significant femtosecond pulse-shaping capabilities were already available in the 
early 1990s, with the development of a programmable multi-element liquid-crystal phase 
modulator that operated on a millisecond time scale [397]. Devices with two LCMs made 
possible simultaneous and independent phase and amplitude modulation of spectral 
components [398,399]. Similar capabilities are also available with AOM-based pulse 
shapers [100]. These and other developments have been reviewed [4,5,12,27]. During 
the last decade, physical and chemical applications of AFC motivated further advances 
in femtosecond pulse-shaping technology, including arbitrary amplitude and phase 
modulation in an acousto-optic programmable dispersive filter [401], enhanced resolution 
of LCMs [402,403], compact and robust pulse-shaping [404], pulse-shape modulation at 
nanosecond time scales using an electro-optical gallium arsenide array with controlled 
waveguides [405], and spectral line-by-line pulse shaping [406,407]. The development 
of polarization pulse shaping [108,409] brought an additional dimension to control of 
quantum phenomena, which is particularly important in some applications (e.g., for 
increasing the yield of multiphoton ionization in molecules); recent improvements in this 
area also include full control of the spectral polarization of ultrashort laser pulses [410], 
simultaneous phase, amplitude, and polarization shaping [411-415], and a simplified 
ultrafast polarization shaper using a birefringent prism [416]. Most recently, the 
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shaping of ultraviolet (UV) femtosecond pulses has been demonstrated, including phase 
modulation [417], simultaneous phase and amplitude modulation [118], and polarization 
shaping [419,420]. 

5.2. Optical applications of AFC 

The AFC approach can be used to produce optical fields with prescribed properties, 
which, in turn, can be applied to control physical and chemical phenomena (e.g., in 
atoms, molecules, and semiconductor structures). In particular, some of the earliest 
AFC experiments aimed at the maximal compression of femtosecond laser pulses 
[421-426]. In these experiments, the light produced through SHG of the shaped pulse 
in a thin nonlinear crystal served as the feedback signal. The SHG yield is directly 
proportional to the intensity of the incident light pulse, and, for pulses with a fixed 
energy, the most intense pulse is the shortest one. AFC-optimized compression of 
broadband laser pulses was also demonstrated using a feedback signal derived from 
two-photon absorption in semiconductors [127]. Application of AFC makes it possible 
to generate maximally compressed laser pulses in a simple and effective way, without 
requiring knowledge of the input pulse's shape. Such adaptive pulse compressors 
(with AOM-based pulse shapers and SHG-based feedback signal) are now employed as 
built-in components in some commercially available femtosecond amplification systems. 
However, since the amplification process is nearly linear, the full-scale application of 
AFC is usually not necessary, as pulses can be compressed in a single feedback step using 
spectral interferometry. The resulting transform-limited pulses can be used as a starting 
point for the study and control of various photophysical and photochemical processes 
(e.g., they can be used to excite and track localized fine-structure and Rydberg wave 
packets in atoms and vibrational wave packets in molecules). In many AFC experiments, 
transform-limited pulses are used as a reference, to separate off the intensity dependence 
which is ubiquitous in nonlinear processes. 

Another optical application of AFC is optimal amplification of chirped femtosecond 
laser pulses [428,429]. The AFC method is used to minimize the higher-order phase 
dispersion that is inherent in the amplification process. Furthermore, AFC was applied 
to generate almost arbitrary target temporal shapes starting with uncharacterized input 
pulses [430,431]. These experiments used a cross-correlation measurement [430] and 
electric field characterization via frequency-resolved optical gating (FROG) [131] of 
the output pulses as the feedback signal. As polarization shaping technology for 
femtosecond laser pulses developed [408,409], AFC was used to generate pulses with 
target polarization profiles [432,433]. One experiment [432] used the SHG feedback 
signal to compensate for material dispersion and time-dependent modulation of the 
polarization state. Another experiment [433] employed a sophisticated feedback signal 
based on dual-channel spectral interferometry to generate shaped femtosecond pulses 
whose ellipticity increased at a constant rate. In a further development, a recent 
experiment [434] used polarization-shaped laser pulses and AFC to manipulate the 
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optical near field on a nanometer scale. 

5.3. AFC of high-harmonic generation 

Among important physical applications of AFC is coherent manipulation of soft X-rays 
produced via high-harmonic generation. In a pioneering experiment, Murnane, Kapteyn, 
and co-workers [435] used shaped ultrashort, intense laser pulses (with 6-8 optical 
cycles) for AFC of high-harmonic generation in atomic gases. Their results demonstrate 
that optimally shaped laser pulses identified by the learning algorithm can improve 
the efficiency of X-ray generation by an order of magnitude, manipulate the spectral 
characteristics of the emitted radiation, and "channel" the interaction between nonlinear 
processes of different orders. All these effects result from complex interferences between 
the quantum amplitudes of the atomic states, created by the external laser field. The 
learning algorithm guides the pulse shaper to tailor the laser field to produce the optimal 
interference pattern. Several consequent AFC experiments [436-438] explored various 
aspects of optimal high-harmonic generation in atomic gases, including the analysis 
of the control mechanism via a comparison of experimental data with predictions of 
theoretical models. Further experimental studies used AFC for optimal spatial control 
of high-harmonic generation in hollow fibers [439,440], optimal control of the brilliance of 
high-harmonic generation in gas jet and capillary setups [441], and optimal control of the 
spectral shape of coherent soft X-rays [ i i'2]. The latter work [442] has been a precursor 
to a more recent development, in which spectrally shaped femtosecond X-ray fields were 
themselves used to adaptively control photofragmentation yields of SFg [443]. Advances 
in optimal control of high-harmonic generation (including related AFC experiments) 
have been recently reviewed [35,41]. 

Beyond the physical interest in achieving control over high-harmonic generation, 
these experiments also demonstrated a dramatic degree of inherent robustness to laser- 
field noise in strongly non-linear control. This behavior can be understood in terms 
of an extensive null space of the Hessian at the top of the control landscape, implying 
a very gentle slope near the global maximum [190, 194]. This characteristic of the 
quantum control landscape makes it possible to tolerate much of the laser noise while 
maintaining a high control yield. Such robustness is expected to be a key attractive 
feature of observable control across virtually all quantum phenomena. 

5.4. AFC of multiphoton transitions in atoms 

Control of bound-to-bound multiphoton transitions in atoms with optimally shaped 
femtosecond laser pulses provides a vivid illustration of the control mechanism based 
on multi-pathway quantum interference. Non-resonant multiphoton transitions involve 
many routes through a continuum of virtual levels. The interference pattern excited 
by the multiple frequency components of the control pulse can enhance or diminish 
the total transition probability. The interference effect depends on the spectral phase 
distribution of the laser pulse. A number of experiments [444-446] used AFC to 
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identify pulse shapes that are optimal for enhancing or cancelling the probability of 
making a transition. In particular, it was possible to tailor dark pulses that do not 
excite the atom at all due to destructive quantum interference. On the other hand, 
AFC was able to find shaped pulses that induce transitions as effectively as transform- 
limited pulses, even though their peak intensities are much lower. Due to the relative 
simplicity of the atomic systems studied, it was possible to compare the results of the 
AFC experiments with theoretical predictions and verify the control mechanism based 
on quantum interference of multiple laser-driven transition amplitudes. In a related 
experiment [447], AFC was helpful for demonstrating that transform-limited pulses are 
not optimal for inducing resonant multiphoton transitions. It was shown that optimally 
shaped pulses enhance resonant multiphoton transitions significantly beyond the level 
achieved by maximizing the pulse's peak intensity. A recent experiment [448] considered 
non- resonant multiphoton absorption in atomic sodium in the strong-field limit. It 
was demonstrated that in this regime the stimulated emission induced by the dynamic 
Stark shift becomes important, which makes transform-limited pulses not optimal for 
strong-field non-resonant multiphoton transitions. AFC was used to discover strong- 
field shaped laser pulses that optimally counteract the dynamic Stark shift-induced 
stimulated emission and thus maximize the absorption probability. 

A more complex problem in atomic physics is control of multiphoton ionization. In 
one experiment [449], AFC was applied to optimize multiphoton ionization of atomic 
calcium by shaped femtosecond laser pulses. The feedback signals were measured using 
ion and electron spectroscopy, and the optimization results were used to elucidate the 
intermediate resonances involved in the photoionization process. Another experiment 
[450] studied photoionization of potassium atoms controlled by phase-locked pairs of 
intense femtosecond laser pulses. Measurements of the Autler-Townes doublet in the 
photoelectron spectra enabled analysis of the induced transient processes. The AFC 
experimental results were helpful for exploring the control mechanism based on the 
selective population of dressed states. 

5.5. AFC of Rydberg wave packets in atoms 

In one of the first applications of AFC, in 1999, Bucksbaum and co-workers [451] 
manipulated the shape of an atomic radial wave function (a so-called Rydberg wave 
packet). Non-stationary Rydberg wave packets were created by irradiating cesium atoms 
with shaped ultrafast laser pulses. A variation of the quantum holography method [452] 
was used to measure the atomic radial wave function generated by the laser pulse. In 
order to reconstruct the wave function, the amplitude of each energy eigenstate in the 
total wave packet was measured independently via state selective field ionization. The 
distance between the measured and target wave packet provided the feedback signal. 
In the weak- field limit, a simple linear relationship exists between the amplitudes of 
the energy eigenstates in the wave packet decomposition and the amplitudes of the 
corresponding spectral components of the control laser field. Based on this relationship. 
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a simple gradient-type algorithm was employed to adjust the spectral phase distribution 
of the control field. AFC equipped with this algorithm was able to change the shape 
of the Rydberg wave packet to match the target within two iterations of the feedback 
control loop. If the wave packet is created in the strong-field regime, then a more 
sophisticated learning algorithm is generally required to implement AFC. 

5. 6. AFC of electronic excitations in molecules 

The first AFC experiment was reported in 1997 by Wilson and co-workers [123]. 
Femtosecond laser pulses shaped by a computer-controlled AOM were used to excite 
an electronic transition in molecules (laser dye IR125 in methanol solution). The 
measured fiuorescence served as the feedback signal in AFC to optimize the population 
transfer from the ground to first excited molecular electronic state. Both excitation 
efficiency (the ratio of the excited state population to the laser energy) and effectiveness 
(the total excited-state population) were optimized. Similar AFC experiments were 
later performed with different molecules in the liquid phase: laser dyes LDS750 in 
acetonitrile and ethanol solutions [453], DCM in methanol solution [454], rhodamine 101 
in methanol solution [455], coumarin 515 in ethanol solution [456], coumarin 6 in a range 
of non-polar solvents (linear and cyclic alkanes) [457], the charge-transfer coordination 
complex [Ru(dpb)3](PF6)2 (where dpb is 4,4'-diphenyl-2,2'-bipyridine) in methanol 
solution [458] and acetonitrile solution [459-461], a donor-acceptor macromolecule (a 
phenylene ethynylene dendrimer tethered to perylene) in dichloromethane solution [462] , 
and perylene in chloroform solution [463]. Two-photon electronic excitations in flavin 
mononucleotide in aqueous solution were controlled using multi-objective optimization 
(a genetic algorithm was employed to simultaneously maximize the fluorescence intensity 
and the ratio of fluorescence and SHG intensities) [464]. Molecular electronic excitations 
were also optimized in AFC experiments in the solid state, with a crystal of a-perylene 
[463,465]. 

Since these AFC experiments are performed in the condensed phase, an important 
issue is the degree of coherence of the controlled dynamics. This issue has been recently 
explored in a series of AFC experiments [457], in which the level of attained control 
was investigated by systematically varying properties of the environment. Specifically, 
AFC was applied to optimize the stimulated emission from coumarin 6 (a laser dye 
molecule) dissolved in cyclohexane, and the recorded optimal pulse shape (characterized 
by a significantly nonlinear negative chirp) was used with several other solvents (linear 
and cyclic alkanes). In these experiments, the molecule was excited in the linear 
absorption regime in order to exclude the trivial intensity dependence characteristic of 
multiphoton processes. The results revealed an inverse correlation between the obtained 
degree of control (as measured by the enhancement of stimulated emission relative to 
that achieved by excitation with the transform-limited pulse) and the viscosity of the 
solvent. This study indicates that the control mechanism involves a coherent process 
(i.e., based on quantum interference of coherent pathways) and that environmentally- 
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induced decoherence limits the leverage of control on the particular molecular system. 
Also, in a recent AFC experiment [462], the photoemission yield of a donor-acceptor 
macromolecule was maximized, and a coherent control mechanism was identified by 
analyzing the pulse optimization process and optimal pulse features. Another AFC 
experiment [466] controlled flavin mononucleotide and riboflavin dissolved in water, 
through a vibronic transition driven by a shaped 400 nm control pulse followed by 
a delayed 800 nm pulse that produced irreversible further excitation. Fluorescence 
depletion served as a feedback signal. The experiment only functioned when the delay 
between the two pulses was less than ~ 1 ps, indicating that the underlying control 
mechanism employs coherent dynamics (see further discussion of this experiment in 
section 5.12). 

5.7. AFC of photodissociation reactions in molecules 

A long-standing goal of photochemistry is selective control of molecular fragmentation. 
During the last decade, AFC employing shaped femtosecond laser pulses has been 
applied to achieve significant successes towards meeting this goal [18,26,37]. Selective 
quantum control of photodissociation reactions in molecules using AFC was first 
demonstrated by Gerber and co-workers in 1998 [124]. They studied photodissociation of 
the organometallic complex CpFe(C0)2Cl (where Cp = C5H5) that contains particular 
types of iron-ligand bonds and exhibits different fragmentation channels upon excitation 
with shaped femtosecond laser pulses. The branching ratio [CpFe(CO)Cl]"'"/[FeCl]"'' was 
maximized and minimized in AFC experiments employing an evolutionary algorithm. 
The experiment was performed in a molecular beam, and the feedback signal was 
obtained from measurements of the ionized photofragments in a time-of-flight mass 
spectrometer. Using AFC, it was possible to change the branching ratio between 5:1 
and 1:1. 

The success of the AFC experiment described above triggered an ongoing wave of 
research activity in this area. In particular, Gerber's group explored various aspects 
of AFC of photodissociation and photoionization reactions in molecules. In one AFC 
experiment [467], the relative yields of photodissociation and photoionization of iron 
pentacarbonyl, Fe(C0)5, were controlled in the gas phase, using femtosecond laser 
pulses with carrier wavelengths at 800 nm and 400 nm (the latter produced via SHG 
of the former). The AFC-based optimization (both maximization and minimization) 
of the branching ratio [Fe(C0)5]"^/Fe"'" demonstrated that the control mechanism is 
not simply intensity-dependent, but rather employs the spectral phase distribution 
of the shaped laser pulse to steer the dynamics of the excited molecular vibrational 
wave packet towards the target reaction channel. Another gas-phase AFC experiment 
[468] also analyzed the relative importance of intensity-dependent and coherent effects 
in control of photochemical reactions that involve nonlinear (multiphoton) optical 
excitations. The control goals were the direct photoionization of CpFe(C0)2Cl (i.e., 
maximization of the [CpFe(C0)2Cl]^ yield) and selective photofragmentation (i.e.. 
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maximization of the branching ratio [CpFe(CO)Cl]"^/[FeCl]"'"). For each pulse shape 
during the AFC-based optimization, the target reaction yield and SHG efficiency 
(which is directly proportional to the pulse intensity) were recorded. In the case 
of direct ionization control, a clear correlation between the [CpFe(C0)2Cl]"'" yield 
and SHG efficiency was observed, which implies that the photoionization control 
mechanism is mainly intensity-dependent. However, for fragmentation control, no 
correlation between the [CpFe(CO)Cl]"'"/[FeCl]"'" ratio and SHG efficiency was found. 
Moreover, for different pulses with the same SHG intensity, a large range of different 
[CpFe(CO)Cl]"'"/[FeCl]"'" values was obtained, depending on the specific pulse shape. 
These results indicate that while photofragmentation involves multiphoton excitation, 
it is not regulated by the pulse intensity alone. Rather, the control mechanism for 
a particular photofragmentation reaction requires a specially tailored laser pulse to 
guide the complex wave packet dynamics towards the desired outcome. Results of 
a similar AFC experiment [469] that analyzed photofragmentation of CH2ClBr in 
the gas phase (including maximization and minimization of the [CH2Br]"'"/[CH2Cl] + 
ratio) also indicate that the control mechanism involves manipulation of the wave 
packet dynamics on neutral dissociative surfaces rather than purely intensity-dependent 
effects. Experiments that demonstrated AFC of photofragmentation in the molecules 
CpFe(C0)2Cl and CpFe(C0)2Br [470] will be discussed later in the context of optimal 
dynamic discrimination of similar quantum systems. 

In 2001, Levis and co-workers [471] used AFC with shaped, strong-field laser pulses 
to demonstrate selective cleavage and rearrangement of chemical bonds in polyatomic 
organic molecules (in the gas phase), including (CH3)2CO (acetone), CH3COCF3 
(trifiuoroacetone) , and CgHsCOCHs (acetophenone). Control over the formation of 
CH3CO from (CH3)2CO, CF3 or CH3 from CH3COCF3, and C6H5CH3 (toluene) from 
C6H5COCH3 was achieved with high selectivity. The use of strong laser fields (with 
intensities of about 10^'^ W/cm^) helps to effectively increase the available bandwidth, 
as transitions to excited molecular states are facilitated by the dynamic Stark shift. 
This effect opens up many reaction pathways which are inaccessible in the weak-field 
(perturbative) regime due to resonant spectral restrictions [18]. While theoretical 
treatment of the complex strong- field molecular dynamics is extremely difficult, this 
complexity in no way affects employment of AFC in the laboratory, where the molecule 
solves its own Schrodinger equation on a femtosecond time scale. By operating at a 
high-duty control cycle, a learning algorithm is typically able to identify optimal laser 
pulses in a matter of minutes. 

Significant attention has been devoted to the analysis of quantum dynamical 
processes involved in molecular photofragmentation control achieved in gas-phase AFC 
experiments with shaped femtosecond laser pulses. Woste and co-workers [172-474] 
studied mechanisms of photofragmentation control for CpMn(C0)3 optimizing the 
branching ratios [CpMn(C0)] + /[CpMn(C0)3]+ and [CpMn(CO)2]+/[CpMn(CO)3] + . 
Weinacht and co-workers [475-478] analyzed mechanisms underlying control of 
photofragmentation in a series of AFC experiments with similar molecules: CH3COCF3 
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(trifluoroacetone), CH3COCCI3 (trichloroacetone) , and CH3COCD3 (tri-deuterated 
acetone). The yield of the [0X3]+ fragment and the [CX3]"''/[CH3]~'" ratio (where 
X is F, CI, and D for trifluoroacetone, trichloroacetone, and tri-deuterated acetone, 
respectively) were optimized in these AFC experiments using intense shaped laser 
pulses. AFC was also used to optimize the branching ratios Br"'"/[CH2Br]+ and 
[CH2l]^/[CH2Br]"'' in photofragmentation of CH2BrI (bromoiodomethane) [478,479]. In 
a number of works [472,474-476], AFC experiments were supplemented by theoretical ab 
initio quantum calculations to help clarify photofragmentation control mechanisms. In 
several other studies [477-479], a change in the basis of the control variables made 
it possible to reduce the dimension of the search space and thus elucidate control 
mechanisms of selective molecular photofragmentation. In another work [480], pump- 
probe spectroscopy was utilized to explore the control mechanism of photofragmentation 
of CHBr2COCF3 (1,1-3,3,3 dibromo-trifiuoroacetone) in AFC experiments with intense 
shaped laser pulses. In particular, optimization of the [CF3]+/[CHBr2]"'" ratio revealed 
a charge-transfer-based control mechanism. 

Selective control of molecular fragmentation in the gas phase was demonstrated 
in several other AFC experiments with shaped femtosecond laser pulses. Woste and 
co-workers [481-484] controlled photoionization and photofragmentation dynamics of 
alkali clusters. Jones and co-workers [185] optimized the S^/S^^ ratios for various 
and M values in strong-field photofragmentation of Sg molecules. It was found that 
optimally shaped pulses dramatically outperform the transform-limited pulses. Wells 
et al. [486] controlled the vibrational population distribution in the transient CO^"*" 
to manipulate the branching ratio of the CO^^ and C"*" -|- O"*" products. Hill and 
co-workers [187] controlled the amplitude of the bending vibrational mode in highly 
ionized CO2 (during strong-field Coulomb explosion) to enhance the symmetric six- 
electron fragmentation channel, COg"*" — )■ 0^^ + C^^ + O^^. They constrained the search 
space by expressing the spectral phase of the laser pulse as a Taylor series, in order to 
elucidate the controlled photodissociation dynamics. Laarmann et al. [488,489] achieved 
selective cleavage of strong backbone bonds in amino acid complexes (in particular, 
a peptide bond in Ac-Phe-NHMe and Ac-Ala-NHMe), while keeping weaker bonds 
intact. Based on these results, they suggested the possibility of employing AFC with 
optimally tailored laser pulses as an analytical tool in mass spectrometry of complex 
polyatomic systems (with potential applicability, e.g., to protein sequencing of large 
biopolymers). In a recent AFC experiment, Levis and co-workers [490] used intense 
laser pulses to manipulate branching ratios of various photofragmentation products 
of dimethyl methylphosphonate (DMMP), a simulant for the nerve agent sarin. The 
optimization in this experiment was performed in the presence of a high background of 
a hydrocarbon and water in the extraction region of a time-of-flight mass spectrometer. 
The ability to achieve highly selective control under these conditions demonstrates that 
AFC may provide the means to identify complex airborne molecules. As mentioned 
in section 5.3, photofragmentation of SFg was controlled (including optimization of 
the ratio [SF5]+/[SF3]+) in an AFC experiment [443] that used spectrally shaped 
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femtosecond X-ray fields produced from intense shaped laser pulses via high-harmonic 
generation. 

Recently, Dantus and co-workers [491, 492] reported a study of molecular 
fragmentation using intense femtosecond pulses, in which they did not use algorithm- 
guided AFC, but rather evaluated a large set of predetermined pulse shapes. They 
concluded that the yields of photofragmentation products are mainly controlled by the 
pulse intensity, while the details of the pulse shape (e.g., the spectral phase distribution) 
are not important. The contradiction between this conclusion and the results of 
numerous AFC experiments discussed above appears to be explained by the fact that 
Dantus and co-workers performed their measurements [191] at laser intensities above 
the saturation threshold for ionization. As pointed out by Levis [19->], the experimental 
conditions employed by the Dantus group obscure connections to coherent processes 
which occur below the saturation threshold. Since these coherent processes play the 
main role in control mechanisms underlying reaction selectivity achieved by optimally 
shaped laser pulses, the proper experimental conditions must be satisfied for coherent 
control of photoinduced molecular fragmentation. Moreover, the search over a set of 
predetermined pulse shapes [491] cannot guarantee discovery of optimal controls, even 
if the employed set is very large. In a parameter space of the size characteristically 
available from a typical pulse shaper, only a dedicated optimization algorithm is capable 
of consistently identifying optimal control fields which commonly lie in a null space of 
the full search space. This situation is fully consistent with classical engineering control 
practice where optimization is the design tool employed almost without exception to 
meet complex operational objectives. 

5.8. AFC of multiphoton ionization in molecules 

The use of polarization-shaped femtosecond laser pulses can significantly enhance the 
level of control over multiphoton ionization in molecules. In 2004, Brixner et al. [194] 
demonstrated that a suitably polarization-shaped laser pulse increased the photoion- 
ization yield in K2 beyond that obtained with an optimally shaped linearly polarized 
laser pulse. This effect is explained by the existence of different multiphoton ionization 
pathways in the molecule involving dipole transitions which are preferably excited by 
different polarization directions of the laser field. Suzuki et al. [495] applied AFC with 
polarization-shaped laser pulses to multiphoton ionization of I2 molecules and optimized 
the production of oddly charged (I^ and and evenly charged (l2^) molecular ions. 
Weber et al. [496] performed AFC experiments with polarization-shaped laser pulses to 
optimize the photoionization yield in NaK molecules. Free optimization of the pulse 
phase, amplitude, and polarization resulted in a higher ionization yield than parame- 
terized optimization with a train of two pulses. 

Woste and co-workers [497-502] investigated AFC of multiphoton ionization in 
K2 and NaK using femtosecond laser pulses with phase and amplitude modulation, 
but without polarization shaping. In particular, good agreement between the optimal 
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shapes of laser pulses obtained theoretically (via QOCT) and experimentally (via AFC) 
was reported [499]. Significant attention was also devoted to studying properties of 
the controlled ionization dynamics and revealing underlying control mechanisms. In 
one series of AFC experiments [500], the photoionization dynamics was explored using 
control pulse cleaning (CPC), which is a process of removing extraneous control field 
features by applying pressure in the optimization algorithm with an appropriate cost 
function on the spectral components of the pulse [378]. Weak pressure applied in the 
case of isotope-selective ionization of K2 was sufficient to remove unnecessary pulse 
components and thus expose the participating vibronic transitions. For ionization of 
NaK, strong pressure was applied and multi-objective optimization was performed. The 
resultant Pareto-optimal curve revealed the correlation of the two conflicting objectives 
of maximizing the ionization yield versus cleaning the control pulse. The optimal 
ionization pathway depends on the CPC strength, which helps to identify the important 
electronic transitions to particular vibrational states. These results demonstrate that 
the spectra of optimal pulses obtained with CPC contain important information about 
the control mechanism. In another series of AFC experiments [501], the control 
mechanism of multiphoton ionization in NaK was analyzed by systematically reducing 
the complexity of the search space. The spectral phase function of the control pulse was 
expressed as a truncated Fourier series, whose parameters were examined with respect 
to the ionization yield and the obtained optimal field. By progressively reducing the 
number of phase modulation parameters, it was possible to generate optimized pulses 
that allowed for a simple mechanistic interpretation of the controlled dynamics. In an 
earlier study, Leone and co-workers [50-3] applied AFC to optimize the weak-field pump- 
probe photoionization signal in Li2 and used first-order time-dependent perturbation 
theory to investigate the dynamics of a rotational wave packet excited by the pump 
pulse and explain the corresponding control mechanism. 

Woste and co-workers [498, 504-506] also demonstrated that AFC is capable of 
achieving isotope-selective ionization of diatomic molecules such as K2 and NaK. They 
showed that optimally tailored control pulses can increase the divergence between the 
dynamics of excited vibrational wave packets in distinct isotopomers (these studies will 
be discussed in more detail in section 5.12 below). 

5.9. AFC of molecular alignment 

The controlled alignment of molecules has attracted considerable attention as it can 
provide a well defined sample for subsequent additional control experiments. At 
high laser intensities (~ 10^^-10^^ W/cm^), dynamical variations of the molecular 
polarization can have a significant effect on alignment. By shaping the temporal 
profile of such an intense femtosecond laser pulse, it is possible to achieve control 
over molecular alignment. Quantum dynamics of laser-induced molecular alignment 
is amenable to theoretical treatment and optimization [216, 507-511], and can be 
successfully controlled using simple ultrafast laser pulses [512-516]. AFC provides a 
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very effective general laboratory tool for alignment manipulation, making the best use 
of the laser resources. AFC of molecular alignment with intense shaped laser pulses was 
successfully demonstrated at room temperature for N2 [517,518] and CO [519]. 

5.10. Applications of AFC in nonlinear molecular spectroscopy 

Shaped femtosecond laser pulses were successfully used to enhance resolution and 
improve detection in several areas of nonlinear spectroscopy and microscopy [11]. Of 
particular interest are experiments that employ AFC to identify optimal pulse shapes. 
One area of nonlinear molecular spectroscopy is control of vibrational modes via 
stimulated Raman scattering (SRS). In the gas phase, a number of AFC experiments 
[520-522] manipulated molecular vibrations excited via SRS by intense ultrafast laser 
pulses in the impulsive regime (i.e., when the duration of the control laser pulse is shorter 
than the vibrational period). In one of these experiments [520], control of the vibrational 
dynamics of K2 was achieved via impulsive SRS in a degenerate four- wave-mixing optical 
setup. Different parameterizations of shaped femtosecond laser pulses in the frequency 
and time domains were employed to decipher the physical mechanism responsible for the 
achieved control. In other gas-phase experiments, mode suppression and enhancement in 
sulfur hexafiuoride [521], mode-selective excitation in carbon dioxide [521], and creation 
of shaped multimode vibrational wave packets with overtone and combination mode 
excitation in CCI4 [522] were demonstrated using impulsive SRS at room temperature 
and high pressures. In the liquid phase, AFC was applied to control relative intensities 
of the peaks in the Raman spectrum, corresponding to the symmetric and antisymmetric 
C-H stretch modes of methanol [523-526]. The modes were excited via SRS in the non- 
impulsive regime (i.e., the duration of the control laser pulse exceeded the vibrational 
period). The control pulse was shaped and the forward scattered Raman spectrum was 
measured to obtain the feedback signal, with the goal of achieving selective control of 
the vibrational modes. However, it was argued [527, 528] that in non-impulsive SRS 
the relative peak heights in the Raman spectrum do not reflect the relative populations 
of the vibrational modes and that control of the spectral features demonstrated in the 
experiments [523-526] does not involve quantum interference of vibrational excitations, 
but rather is based on classical nonlinear optical effects. 

Another important area of nonlinear molecular spectroscopy is control of molecular 
vibrational modes via CARS. Materny and co-workers [529-534], Zhang et al. [535] and 
von Vacano et al. [536] used AFC in a CARS setup to optimally control vibrational 
dynamics in complex molecules. The Stokes pulse was shaped and the feedback 
signal was derived from the intensities observed in the CARS spectrum. In an AFC 
experiment with polymers (fiuorobenzene sulfonate diacetylenes) , selective excitation of 
one vibrational ground-state mode and suppression of all other modes was achieved, 
and the decay times of different modes were modified [529]. In liquid-phase AFC 
experiments, selective enhancement or suppression of one or more vibrational modes 
was demonstrated for toluene [530, 531], benzene [535], and /3-carotene in hexane 
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solution [532]. In a related experiment [533], AFC was used to obtain molecule- 
specific CARS spectra from a mixture of benzene and chloroform. Molecule-specific 
enhancement or suppression of the CARS spectral lines achieved in this AFC experiment 
is an example of optimal dynamic discrimination that will be discussed in more detail 
in section 5.12 below. Another related study [534] compared selective excitation of 
molecular vibrational modes (achieved in AFC experiments with shaped femtosecond 
pulses) in the gas phase with carbon disulfide and the liquid phase with toluene. 
Interestingly, it was found that the relative intensities of the CARS spectral lines 
could be changed more effectively in the liquid phase than in the gas phase. An 
experiment reported by von Vacano et al. [536] employed the method of single-beam 
CARS spectroscopy with shaping of broadband pulses from a photonic crystal fiber. 
Each broadband pulse provides numerous pairs of pump and Stokes frequencies, and the 
spectral phase of the pulse was optimized with AFC to produce the desired interference 
pattern of the molecular vibrational modes. The optimally compressed and shaped 
pulses enabled the unambiguous assignment of the participating vibrational modes of 
toluene between 500 and 1000 cm~^ in a blue-shifted CARS signal. 

Control of molecular vibrational dynamics is possible not only through Raman-type 
processes, but also directly in the infra-red (IR) regime. Zanni and co-workers [537] 
demonstrated selective control of vibrational excitations on the ground electronic state 
of W(C0)6, using shaped femtosecond mid-IR (5.2 /xm, 1923 cm~^) pulses. The spectral 
phase distribution of the pulse was optimized using AFC to achieve selective population 
of the excited vibrational levels of the Tiu CO-stretching mode. Systematic truncation 
of optimal pulses was employed to analyze the control mechanism. In a related AFC 
experiment [538], polarization-shaped mid-IR pulses were used to selectively control 
vibrational excitations of the two carbonyl stretching modes in Mn(C0)5Br. 

5.11. Applications of AFC in multiphoton microscopy 

An important application of AFC with shaped femtosecond laser pulses is in multiphoton 
excited fluorescence (MPEF) microscopy. For a given pulse energy, the transform- 
limited pulse has the maximum peak intensity, which helps to increase the fluorescence 
signal intensity, but unfortunately also increases the rate of photobleaching of the 
molecules (which is especially undesirable with samples of live cells). Use of optimally 
shaped pulses instead of a transform-limited pulse can reduce the bleaching rate, 
enhance spatial resolution, and increase contrast in biological fluorescence imaging. 
In a series of AFC experiments with shaped laser pulses, Midorikawa and co-workers 
[539-542] optimally controlled MPEF microscopy in different fluorescent biomolecules. 
Attenuation of photobleaching by a factor of four (without decreasing the fluorescence 
signal intensity) was demonstrated in two-photon excitation fluorescence (TPEF) from 
a green fluorescent protein [539]. Another AFC experiment [540] achieved selective 
control of two-photon and three-photon fluorescence in a mixture of two biosamples. 
The use of optimally tailored pulses helped to minimize the harmful three-photon 
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fluorescence from the amino acid L-Tryptoplian, witliout a significant loss of useful 
two-photon fluorescence from a green fluorescence protein. Optimally shaped super- 
continuum pulses from a microstructure fiber were used in TPEF microscopy in another 
experimental study [511]. The pulse was shaped prior to propagation through the fiber, 
and AFC maximized the fluorescence signal contrast between two fluorescent proteins. 
A novel phase modulation technique for ultra-broadband laser pulses was developed for 
selective excitation of multiple fluorophores in TPEF microscopy [542]. This technique 
was applied to dual-color imaging of cells containing two types of fluorescent proteins, 
and AFC was employed to find the phase modulation that maximizes or minimizes the 
individual TPEF intensity from one of the fluorophores. 

5.12. Applications of AFC for optimal dynamic discrimination 

Discrimination of similar systems is important for many practical problems in science 
and engineering. In particular, selective identification of target molecules in a mixture 
of structurally and spectroscopically similar compounds is a challenge in such areas 
as selective excitation of multiple fluorescent proteins in microscopy of live samples, 
targeted component excitation in solid-state arrays, and selective transformation of 
chemically similar molecules. Theoretical studies [219,543-546] indicate that quantum 
systems differing even very slightly in structure may be distinguished by means of their 
dynamics when acted upon by a suitably tailored ultrafast control field. Such optimal 
dynamic discrimination (ODD) can in principle achieve dramatic levels of control, and 
hence also provides a valuable test of the fundamental selectivity limits of quantum 
control despite noise and constrained laser resources. AFC provides a very effective 
laboratory means for practical implementation of ODD. 

In 2001, Gerber and co-workers [547] experimentally demonstrated selective multi- 
photon excitation of two complex molecules, a laser dye DCM and [Ru(dpb)3](PF6)2, 
in methanol solution. The goal was to electronically excite DCM while simultaneously 
suppressing electronic excitation of [Ru(dpb)3]^"'". While these two molecules are elec- 
tronically and structurally distinct, the DCM/[Ru(dpb)3]^"'" emission ratio is practically 
unaffected by variations in single control parameters such as wavelength, intensity, and 
linear chirp. Nevertheless, selective excitation was successfully achieved using AFC with 
shaped femtosecond laser pulses. The DCM/[Ru(dpb)3]^"'' emission ratio was used as 
the feedback signal, and the evolutionary algorithm identified optimally shaped con- 
trol pulses that improved the signal by approximately 50%. These results obtained in 
the presence of complex solvent/solute interactions aroused significant interest to ODD. 
As mentioned in section 5.7, Gerber and co-workers [470] demonstrated AFC of pho- 
toproduct branching ratios in CpFe(C0)2Cl and CpFe(C0)2Br. Despite the chemical 
similarity of these two molecules, AFC was sensitive enough to detect differences due to 
the electronic metal-halogen bonding properties. This finding suggests the possibility of 
performing ODD of individual compounds in mixtures of chemically similar molecules. 
Other examples of ODD include molecule-specific manipulation of CARS spectra from 
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a mixture of benzene and chloroform [533], selective excitation of multiple fiuorophores 
in TPEF microscopy [541,542], and quantitative differentiation of dyes with overlapping 
one-photon spectra [548]. 

A recent experimental demonstration of ODD by Roth et al. [KKi] achieved dis- 
tinguishing excitations of two nearly identical flavin molecules in aqueous phase. The 
absorption spectra for flavin mononucleotide (FMN) and riboflavin (RBF) are practi- 
cally indistinguishable throughout the entire visible and far UV. This implementation of 
ODD used a shaped UV pulse centered at 400 nm and a time-delayed unshaped IR pulse 
centered at 800 nm. The first pulse creates a coherent vibrational wave packet on an ex- 
cited electronic state, and the second pulse disrupts the wave packet motion and results 
in additional excitation to a higher state and consequential depletion of the recorded 
fluorescence signal. The effect of slight differences in the vibronic structure of the two 
molecules upon the dynamics of the excited wave packets is amplified by tailoring the 
spectral phase of the UV pulse. Since further excitation produced by the second pulse 
depends on the precise structure, position, and coherence of the tailored wave packet, 
it is possible to dynamically interrogate the two statically nearly identical systems and 
thereby produce a discriminating difference in their respective depleted fluorescence sig- 
nals. In contrast, if the UV pulse is transform- limited, then the fluorescence depletion 
signals from the flavins are indistinguishable. UV pulse shapes that optimally discrimi- 
nate between FMN and RBF were discovered using AFC. The optimized depletion ratio 
-Dfmn/-Drbf could be changed by ~ ±28%, despite the initially indistinguishable lin- 
ear and nonlinear optical spectra. Although the laser resources consisted of a modest 
~ 3.5 nm of UV bandwidth and ~ 10 nm of IR bandwidth, significant selectivity was 
achieved with optimal UV pulses working in concert with the time-delayed unshaped 
IR pulse. Although the static spectra appear nearly identical, subtle differences in 
the molecular structure are nonetheless made profound during the tailored evolution of 
wave packets generated by optimal controls. System complexity (e.g., high vibrational 
state density, thermal population, solvent-induced line broadening) effectively enhances 
the ODD capabilities of the control field and compensates for the limited bandwidth 
constraint, thus making dramatic levels of control possible even in the weak-field limit. 

Another example of ODD is isotope-selective ionization of molecules demonstrated 
in a number of AFC experiments by Woste and co-workers [498,504-506]. In particular, 
in an illustrative study [505,506], they applied shaped femtosecond laser pulses to the 
39,39j<^^ and '^^'^^K2 isotopomers and used AFC to maximize and minimize the isotope ion 
ratio R = /(^^'^^K2)//(^^'^^K2). K2 molecules can be ionized in a three-photon process 
at relatively low pulse energies within the available wavelength range (810-833 nm). 
Differences between the evolutions of vibrational wave packets on an excited electronic 
state in the two isotopomers can be amplified by optimized control fields. Operation in 
this fashion made it possible to achieve significant selectivity of isotope ionization, with 
a variation by a factor of -RmaxZ-Rmin ^ 140 between the maximal and minimal values 
of the isotope ion ratio. Information about the dynamics of the controlled vibrational 
wave packets was extracted from the optimal pulse shapes to help reveal ionization 
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pathways [500,505,506]. 

5.13. AFC of energy flow in biomolecular complexes 

Applications of quantum control to increasingly complex molecular systems have 
been considered. In particular, Motzkus and co-workers used AFC with shaped 
femtosecond laser pulses to control and analyze the energy flow pathways in the light- 
harvesting antenna complex LH2 of Rhodopseudomonas acidophila (a photosynthetic 
purple bacterium) [549,550] and /9-carotene [551]. They demonstrated that by shaping 
the spectral phase distribution of the control pulse, it is possible to manipulate the 
branching ratio of energy transfer between intra- and inter-molecular channels in the 
donor-acceptor system of the LH2 complex [549]. Analysis of the transient absorption 
data was used to decipher the control mechanism and identify the molecular states 
participating in energy transfer within LH2 [550] and /3-carotene [551]. These results 
indicate that coherent quantum control is possible even in very complex molecular 
systems in a condensed-phase environment. 

5.14. AFC of photoinduced electron transfer 

AFC has been applied to quantum control of inter-molecular photoinduced electron 
transfer. Yartsev and co-workers [552] reported an AFC experiment that maximized 
the yield of ultrafast electron injection from the sensitizer to Ti02 nanocrystals in the 
core part of a dye-sensitized solar cell. The electron transfer process was monitored 
using the transient absorption signal. The impulsive structure of the optimal laser pulse 
was observed to correlate with the coherent nuclear motion of the photoexcited dye. 
The pulse shape and the transient absorption kinetics were explained by an impulsive 
stimulated (anti-Stokes) Raman scattering process, followed by electronic excitation. 

5.15. AFC of photoisomerization in complex molecules 

The control of molecular structure transformations is a coveted goal in chemistry. In 
particular, control of cis-trans isomerization has attracted much attention due to the 
importance of this process in chemistry and biology (e.g., it is a primary step of 
vision). AFC of cis-trans photoisomerization in cyanines (in the liquid phase) with 
shaped femtosecond laser pulses was first reported by Gerber and co-workers [553]. This 
experiment demonstrated that by using optimally shaped laser pulses it is possible to 
enhance or reduce isomerization efficiencies. The mechanism underlying isomerization 
control in this experiment was discussed in a number of theoretical works [554-556]. 
In particular, Hoki and Brumer [551] suggested that isomerization control involves 
an incoherent pump-dump process and that the role of quantum coherence effects 
in the evolution of the excited vibrational wave packet is negligible due to strong 
environmentally-induced decoherence. On the other hand. Hunt and Robb [555] and 
Improta and Santoro [556] used a more sophisticated model and argued that control of 
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isomerization does rely on quantum coherence of the photoexcited vibrational wave 
packet that moves on the multidimensional potential-energy surface of an excited 
electronic state. Moreover, Hunt and Robb [555] showed that the generation of the 
trans versus cis product is affected by the presence of an extended conical intersection 
seam on the potential-energy surface, and that photoisomerization can therefore be 
coherently controlled by tuning the distribution of momentum components in the 
photoexcited vibrational wave packet. The validity of this coherent control mechanism 
was corroborated in a further AFC experiment by Yartsev and co-workers [557]. They 
demonstrated that optimally shaped laser pulses can be used to modify the momentum 
composition of the photoexcited wave packet and thus achieve significant control of the 
absolute yield of isomerization (i.e., the photoisomer concentration versus the laser pulse 
energy). The coherent character of liquid-phase control of cis-trans photoisomerization 
in cyanines was further studied in another AFC experiment by Yartsev and co-workers 
[558]. They used a control scheme with an unshaped pump pulse and a time-delayed 
shaped dump pulse (an unshaped probe pulse was also applied to measure the effect 
of control). By using the optimally shaped dump pulse, they achieved control of 
photoisomerization closer to the decisive points of the reaction. This approach made it 
possible to explore the effect of the wave-packet's momentum composition at different 
time scales and assign the dynamics to distinct parts of the excited-state potential. 

AFC of the retinal molecule in bacteriorhodopsin (from the all-trans to the 13- cis 
state) was demonstrated by Miller and co-workers [559] . This experiment employed both 
phase and amplitude modulation of femtosecond laser pulses and operated in the weak- 
field regime (with pulse energies of 16-17 nJ). By using optimally shaped pulses, it was 
possible to enhance or suppress the quantity of molecules in the 13- cis state by about 
20%, relative to the yield observed using a transform- limited pulse with the same energy. 
They further explored the mechanism of coherent control of retinal photoisomerization in 
bacteriorhodopsin using time- and frequency-resolved pump-probe measurements [560]. 
Experimental data together with a theoretical analysis suggest that the isomerization 
yield depends on the coherent evolution of the photoexcited vibrational wave packet 
on an excited-state potential-energy surface in the presence of a conical intersection. 
According to this analysis, control of retinal photoisomerization is dominated by 
amplitude modulation of the spectral components of the excitation pulse. Gerber and 
co-workers [5(1 1] also demonstrated AFC of retinal isomerization in bacteriorhodopsin. 
In this experiment, they pioneered the control scheme with unshaped-pump and time- 
delayed shaped-dump femtosecond laser pulses. As mentioned above, the use of the 
optimally shaped dump pulse makes it possible to control the molecule in a region 
of the potential-energy surface where the decisive reaction step occurs. Moreover, by 
changing the time delay between the pulses, it is possible to obtain information on the 
wave packet evolution. 

The role of quantum coherence effects in control of retinal isomerization in 
bacteriorhodopsin is still not fully clear, as a recent experiment by Bucksbaum and co- 
workers [562] found no dependence of the isomerization yield on the control pulse shape 
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at pulse energies below 30 nJ. In the high-intensity regime (at pulse energies above 30 
nJ), they found that the yield of the 13- cis isomer is maximized by a transform-limited 
pulse, which could indicate that the yield depends only on the pulse intensity, with 
quantum coherence not playing a significant role. These findings (especially for lower 
intensities) apparently contradict the optimization results obtained by Miller and co- 
workers [559]. It is possible that these discrepancies could be explained by differences 
in experimental setups. In particular, Bucksbaum and co-workers [562] used only phase 
modulation of the control pulse, whereas Miller and co-workers [559,560] argued that 
control is mainly achieved by amplitude modulation. Additional experimental and 
theoretical work will be needed to fully explore the mechanisms underlying condensed- 
phase control of photoisomerization in complex molecular systems and clarify the role of 
quantum coherence in the controlled dynamics. For example, recent theoretical studies 
[563, 564] suggest that coherent control of photoisomerization and other branching 
reactions in an excited state may be affected by and, moreover, take advantage 
of environmentally-induced relaxation effects. Such cooperation between coherent 
control and environmentally-induced decoherence may be important in various quantum 
phenomena [565] and hence its optimal exploitation deserve further investigation. 

Other examples of structural transformations in complex molecules include ring 
opening in cyclohexadiene along with isomerization as well as cyclization reactions in cis- 
stilbene. Carroll et al. [566, 567] demonstrated AFC of the photoinduced ring-opening 
reaction of 1,3- cyclohexadiene (CHD) to form 1,3,5-cis-hexatriene (Z-HT). The feedback 
signal was obtained from measurements of the UV absorption spectrum. The learning 
algorithm was able to identify optimal pulse shapes that increased the formation of Z- 
HT by a factor of two. For a different control objective, the AFC optimization produced 
pulse shapes that decreased solvent fragmentation while leaving the formation of Z-HT 
essentially unaffected. Kotur et al. [568] used AFC with shaped ultrafast laser pulses 
in the deep UV to control the ring opening reaction of CHD to form 1,3,5-hexatriene 
(HT). The feedback signal was obtained from measurements of fragmentation products 
following strong-field ionization with a time-delayed IR laser pulse. The learning 
algorithm discovered shaped UV pulses that increased the HT yield by ~ 37% relative to 
an unshaped (nearly transform-limited) pulse of the same energy. Greenfield et al. [569] 
demonstrated AFC of the photoisomerization and cyclization reactions in cis-stilbene 
dissolved in n-hexane. This experiment employed phase-modulated 266 nm femtosecond 
pulses to maximize or minimize the yields of cis- to trans-stilbene isomerization as well as 
cis-stilbene to 4a,4{)-dihydrophenanthrene cyclization. The yields of both isomerization 
and cyclization were minimized by transform-limited pulses that enhanced competing 
multiphoton processes, while the yields were maximized by complex pulse shapes that 
helped to avoid multiphoton effects. 



Control of quantum phenomena: Past, present, and future 



42 



5.16. AFC of nuclear motion in fullerenes 

FuUerenes are a class of molecules of considerable interest in many areas of science. 
Laarmann et al. [570] employed AFC-optimized femtosecond laser pulses to coherently 
excite large-amplitude oscillations in Cgo fullerene. The structure of the optimal pulses 
in combination with complementary two-color pump-probe data and time-dependent 
density functional theory calculations provided information on the underlying control 
mechanism. It was found that the strong laser field excites many electrons in Ceo, and 
the nuclear motion is excited, in turn, due to coupling of the electron cloud to a radially 
symmetric breathing mode. Despite the complexity of this multi-particle system with 
various electronic and nuclear degrees of freedom, the optimal control fields generated 
essentially one-dimensional oscillatory motion for up to six cycles with an amplitude of 
~ 130% of the molecular diameter. 

5.17. Applications of AFC in semiconductors 

Quantum control has found applications beyond atomic and molecular phenomena. 
In particular, it is possible to use optimal control methods to manipulate various 
processes in semiconductors. Kunde et al. [571,572] demonstrated AFC of semiconductor 
nonlinearities using phase- modulated femtosecond laser pulses, with the purpose of 
creating an ultrafast all-optical switch. The feedback signal was obtained by measuring 
the differential transmission (DT) in a control-probe setup. Optimizations were 
performed on the spectrally integrated DT as well as DT in narrow spectral windows. 
The learning algorithm was able to identify optimal pulse shapes that enhanced ultrafast 
semiconductor nonlinearities by almost a factor of four. Chung and Weiner [573] 
used AFC with phase-modulated femtosecond laser pulses to coherently control two- 
photon-induced photocurrents in two different semiconductor diodes. Because of their 
spectrally distinct two-photon absorption responses, the diodes generated noticeably 
different photocurrent yields depending on the pulse shape. An evolutionary algorithm 
guided the AFC experiment to discover pulse shapes that maximize or minimize the 
photocurrent yield ratio for the two diodes. 

5.18. AFC of decoherence 

Manipulation of quantum interference requires that the system under control remains 
coherent, avoiding (or at least postponing) the randomization induced by coupling 
to an uncontrolled environment. Therefore, the ability to manage environmentally- 
induced decoherence would bring substantial advantages to control of many physical and 
chemical phenomena. In particular, decoherence is a fundamental obstacle to quantum 
information processing [175], and therefore the ability to protect quantum information 
systems against decoherence is indispensable. 

The possibility of using AFC for optimal suppression of decoherence was first 
proposed by Brif et al. [571], and numerical simulations in a model system were 
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performed by Zhu and Rabitz [387]. Walmsley and co-workers [118] have recently 
used AFC to achieve coherent control of decoherence of molecular vibrational wave 
packets in the laboratory. The concept underlying this experiment is the use of coherent 
preparation of the quantum system to alter non-unitary decoherent dynamics induced 
by an uncontrolled environment. In this experiment, a gas-phase ensemble of K2 at 
400 °C was irradiated by a shaped femtosecond laser pulse, inducing a vibrational wave 
packet in the lowest excited electronic state of the molecules. This wave packet 

undergoes dephasing (a form of decoherence that does not involve dissipation of energy). 
Dephasing is caused by coupling of the vibrational mode to the thermalized rotational 
quasi-bath. The amplitude of quantum beats in the fluorescence signal (measured at 
a chosen delay time after the excitation pulse) served as the feedback signal. This 
amplitude provides an estimate of the degree of wave packet localization in the phase 
space and therefore can be used as a coherence surrogate. The optimal pulse identified 
by AFC increased the quantum-beat visibility from zero to more than four times the 
noise level and prolonged the coherence lifetime by a factor of ~ 2 relative to the 
beats produced by the transform-limited pulse. The main characteristic of the optimal 
pulse is a high degree of linear negative chirp. This indicates that the mechanism of 
decoherence control is based on exciting a wave packet which is initially amplitude- 
squeezed in the phase-space representation. A theoretical analysis confirmed that the 
coherence hfetime is extended if the initial state is amplitude-squeezed (corresponding 
to the initial orientation in the phase space along the classical trajectory). 

A well-known strategy for suppressing decoherence in quantum systems is through 
application of pulse sequences designed to dynamically decouple the system from the 
environment [575-587]. Experiments [588-593] have employed theoretically designed 
pulse sequences based on particular environment models (i.e., an example of open- 
loop control). However, the actual noise power spectra for realistic environments can 
significantly differ from the models. To overcome this difficulty, Bollinger and co- 
workers [594] recently used AFC to tailor the dynamical decoupling pulse sequence 
to an actual experimental noise environment. In this experiment, the system was a 
quantum memory realized in an array of ~ 1000 ^Be"*" ions (cooled to a temperature 
of ~ 1 mK) in a Penning ion trap, with qubit states realized using a ground-state 
electron-spin-flip transition. These qubit states are highly susceptible to magnetic fleld 
fluctuations, making such noise a signiflcant source of decoherence. The qubits were 
coherently controlled by directly driving the ~ 124 GHz transition in a microwave setup 
similar to optical ones. A sequence of microwave tt pulses used for qubit control in 
this laboratory conflguration is technologically quite different from shaped femtosecond 
optical laser pulses typically employed in molecular control experiments; however, the 
fundamental concept of AFC is still fully applicable. The feedback signal was obtained 
from fluorescence detection on a cycling transition (with decoherence-induced errors 
manifested as non-zero fluorescence). The Nelder-Mead simplex method was utilized 
to search for optimal pulse positions in a flxed-length sequence of n pulses. Optimal 
pulse sequences discovered in the AFC experiment, without a priori knowledge of the 
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noise environment, suppressed the qubit error rate by a factor of five to ten relative to 
benchmark model-based sequences. 

5.19. Algorithmic advances for laboratory AFC 

The learning algorithm is an important component of laboratory AFC. The majority of 
AFC experiments employ stochastic algorithms such as evolutionary strategies [595] and 
genetic algorithms [590]. Historically, genetic algorithms were characterized by the use 
of recombination, while evolutionary strategies primarily relied on mutation; however, 
modern algorithms guiding AFC experiments and simulations typically incorporate both 
types of genetic operations and are variably called genetic algorithms or evolutionary 
algorithms. These algorithms are very well suited to laboratory optimizations as they 
naturally match the discrete structure of control "knobs" (e.g., the pixels of a pulse 
shaper) and are robust to noise. Moreover, robustness to noise in AFC experiments 
can be enhanced by incorporating the signal-to-noise ratio into the control objective 
functional [-378,597]. Various aspects of evolutionary algorithms and their application 
to AFC of quantum phenomena were assessed [598,599]. Evolutionary algorithms can be 
also used in multi-objective optimization [600,601], and the application of this approach 
to quantum control problems was studied theoretically [216-218,602] and demonstrated 
in AFC experiments [464,500]. 

Other types of stochastic algorithms include, for example, simulated annealing [603] 
and ant colony optimization [604, 605]. Simulated annealing was utilized in some 
AFC experiments [430,542], and it seems best suited to situations where just a few 
experimental parameters are optimized [606,607]. Ant colony optimization recently has 
been used in an AFC simulation [390], but it is yet to be tested in the laboratory. 

As mentioned in section 4, the absence of local traps in control landscapes for 
controllable quantum systems has important practical consequences for the optimization 
complexity of AFC experiments. In particular, deterministic search algorithms can be 
used to reach a globally optimal solution. Deterministic algorithms (in particular, the 
downhill simplex method) were successfully implemented in several AFC experiments 
[431,569,594]. Recently, Roslund and Rabitz [371] demonstrated the efficiency of a 
gradient algorithm in laboratory AFC of quantum phenomena. They implemented 
a robust statistical method for obtaining the gradient on a general quantum control 
landscape in the presence of noise. The experimentally measured gradient was 
utilized to climb along steepest-ascent trajectories on the landscapes of three quantum 
control problems: spectrally filtered SHG, integrated SHG, and excitation of atomic 
rubidium. The optimization with the gradient algorithm was very efficient, as it 
required approximately three times fewer experiments than needed by a standard genetic 
algorithm in these cases. High algorithmic efficiency is especially important for AFC of 
laser-driven processes in live biological samples, as damage (e.g., due to photobleaching) 
can be reduced by decreasing the number of trial laser pulses. Still, evolutionary or other 
stochastic algorithms may be preferable over deterministic algorithms in many AFC 
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experiments due to their inherent robustness to noise. Hybrid stochastic-deterministic 
algorithms (e.g., derandomized evolution strategies) seem to offer the most ffexibility 
and efficiency [608]. 

6. The role of theoretical quantum control designs in the laboratory 

A very significant portion of theoretical research in the area of quantum control is 
devoted to model-based computations which employ QOCT (or other similar methods) 
to design optimal control fields for various physical and chemical problems. Such 
computations are widespread in theoretical studies of molecular applications of quantum 
control and are becoming increasingly popular in considering control of quantum 
information systems, including optimal operation of quantum gates and optimal 
generation of entanglement (see section 3.6). Notwithstanding these extensive QOCT- 
based control field designs, experiments seeking optimal control of molecular processes 
overwhelmingly employ AFC methods as described in section 5. Such experiments in 
most cases work remarkably well with random initial trials, and thus exhibit no evident 
need to operate or possibly start with theoretical control designs. This raises important 
questions about the practical usefulness of open-loop control and role of theoretical 
methods such as QOCT in control experiments [23]. In considering this matter it is 
important to keep in mind that the AFC procedure grew out of observations from 
QOCT simulations and associated analyses. 

6.1. The effect of system complexity 

As discussed in section 3.7, the practical laboratory relevance of theoretical designs 
depends on the complexity of the system of interest, with simpler cases yielding 
theoretical models closer to reality. For example, in numerous NMR experiments 
employing RF fields to manipulate nuclear spins [153-157], including NMR realizations 
of quantum gates [159-162], theoretically designed sequences of pulses (some of which 
were developed using QOCT [163, 164]) can function quite well. The model of a 
collection of spins (with empirical coupling and decay constants) interacting with 
classical fields is often sufficiently accurate for NMR-based applications, allowing for 
successful employment of open-loop control. A QOCT-based design was also successfully 
applied experimentally to enhance robustness of single-qubit gate operations in a system 
of trapped ions [353]. 

At the other extreme of system complexity are electronic and vibrational processes 
in polyatomic molecules whose dynamics cannot be accurately modelled at the present 
time. An objective assessment is that models used for polyatomic molecules in control 
computations are currently too simplified and computational techniques inadequate 
for the true levels of complexity, resulting in theoretical designs that are not directly 
applicable to control experiments which work with real systems. There are several 
aspects of laser control of molecules, which make the difference between theoretical 
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models and actual systems very important. First, optimal control is generally based 
on creating interference of multiple quantum pathways, which can be very sensitive 
to the detailed properties of the system (e.g., evolution of a laser- induced vibrational 
wave packet can be strongly influenced by small variations of a molecular potential 
energy surface, as evident from experimental ODD results [466]). Therefore, even small 
inaccuracies of theoretical models or associated numerical procedures may result in 
control designs that are not optimal for the actual systems. Second, the wavelength 
of the laser field is typically much larger than the molecule, which makes it impossible 
to focus the field on a particular group of atoms. In many applications, the goal is to 
attain control (e.g., generate an excitation or break a bond) in a localized part of a 
molecule. However, the external control field interacts with the molecule as a whole, 
possibly exciting a multitude of accessible transitions within the bandwidth of the laser 
(including multiphoton processes). The "global" effect of photonic reagents becomes 
ever more important with stronger fields for which resonance-based "localization" is not 
valid. An ab initio quantitative theoretical account of laser-driven molecular dynamics 
is presently not feasible, unless the studies are limited to cases of very simple molecules 
and weak fields. Third, it is difficult to calibrate the laser and pulse-shaping apparatus to 
reliably reproduce theoretical control designs in the laboratory. In many cases, directly 
using AFC optimization is much easier and much more effective, than calibrating the 
laser and pulse-shaper for generation of theoretically computed control fields with the 
required accuracy. 

These considerations lead to the conclusion that open-loop control experiments 
employing theoretical designs may be useful for some systems and impractical for others, 
depending on how well the system is known and which computational capabilities are 
available, consequently determining how accurately the controlled dynamics can be 
modelled. Thus, the boundary between the systems for which modelling is sufficiently 
reliable and the systems for which it is not, depends on available Hamiltonian data, 
numerical algorithms and computational power. Of course, with time, better modelling 
will become available for more complex systems, although the exponential increase of 
the Hilbert-space dimension with the system complexity is a fundamental feature of 
quantum mechanics, which significantly hinders the effectiveness of numerical control 
designs for practical laboratory implementation. 

Consider, for example, quantum information processing systems which are typically 
modelled as collections of two-level particles (qubits) with controlled interactions 
between them [175]. Prima facie, such a system appears to be quite simple, so 
that controls for all desired transformations can be theoretically designed (e.g., using 
QOCT). However, the difficulty of accurately modelling the actual environmental 
noise is significant even for simple few-qubit systems. As was recently demonstrated 
with trapped-ion qubits, dynamical decoupling pulse sequences obtained via AFC 
significantly outperformed the best available theoretical designs [591]. Moreover, as 
the Hilbert-space dimension increases exponentially with the number of qubits, the 
unwanted effects of uncontrolled couphngs between the qubits in multi-particle systems 
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will be very difficult to model and, consequently, to manage via open-loop control. 
Therefore, the effectiveness of theoretical control designs for realistic quantum computers 
will depend on the ability to engineer systems in which couplings between small blocks 
of qubits can be made arbitrarily small. Hopefully, further technological advances will 
make such systems available, which in turn will make open-loop control with QOCT- 
based designs useful for practical quantum computing. 

6.2. Importance of theoretical control designs for feasibility analysis 

In molecules, the interactions between the atoms are inherently strong in order to hold 
the atoms together. Therefore, in the foreseeable future, for optimal manipulation 
of electronic and vibrational processes in molecules with four and more atoms, AFC 
will continue to be much more effective than employing theoretical control designs. 
Notwithstanding this assessment, theoretical control studies should continue to have 
high significance; however, for complex systems the value of theoretical studies is not 
in generating specific control designs for immediate laboratory use. Control solutions 
obtained via theoretical model-based computations (in particular, those employing 
QOCT) should play an important role by advancing the general understanding of the 
character of controlled dynamics and control mechanisms. One practically important 
issue is that while each cycle of a typical AFC experiment is very fast (from microseconds 
to milliseconds) and cheap, the initial setup of the experiment is usually quite difficult 
and expensive, since advanced methods of pulse shaping and control-yield measurement 
need to be incorporated together as well as adjusted to the particular nature of a physical 
or chemical system. Therefore, theory can be especially important in exploring the 
feasibility of various control outcomes for quantum dynamics of a system of interest; even 
semiquantitative modelling may be successful for such purposes in many applications. 
Theoretical control simulations can provide important guidance for the selection of the 
experimental configuration, thereby helping to make AFC a more effective practical 
tool. Additional such high value utilizations of theory and simulations can be expected 
in the future. 

6.3. Open-loop quantum control experiments with non-optimal designs 

The open-loop control procedure is not limited to the use of optimal theoretical designs 
generated via QOCT and similar methods. Moreover, optimality in not always required 
in quantum control. In the conceptually allied field of synthetic chemistry, progress has 
often been achieved via intuition-guided trials leading to a gradual increase of reaction 
yields. Following this venerable tradition, some recent open-loop control experiments 
seek improvement by employing ultrafast shaped laser pulses with so-called "rational" 
or "judicious" control designs obtained using a combination of intuition and arguments 
based on some knowledge of system properties (e.g., spectral information or symmetry). 
This approach is popular in nonlinear spectroscopy and microscopy [44,609-622] as well 
as some other atomic [623-631], molecular [512-516,632,633], and solid-state [634-637] 
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applications. While in many situations such "rationally" designed control pulses enhance 
the spectroscopic resolution or increase the control yield as compared to results obtained 
with transform-limited pulses, in general they are not optimal. Experience gained 
from numerous quantum control experiments indicates that intuition generally fails to 
discover the most effective controls (except for the simplest systems), and therefore in 
most cases the degree of control can be increased via closed-loop optimization employing 
a suitable learning algorithm. In some situations, intuition-driven control may be 
effective for providing a guide to initial fields for subsequent optimization under AFC. 

7. Concepts and applications of real-time feedback control 

Feedback is very important in classical engineering where it is routinely used for control 
of complex systems in the presence of uncertainties. In quantum control, two important 
approaches based on the concept of feedback have been introduced for similar reasons. 
One is AFC, which was extensively discussed in section 5. A fundamental characteristic 
of AFC is that in each control cycle a fresh quantum ensemble is used (either a new 
sample is prepared or the system is reset to its initial state before each run), which 
makes measurement back action irrelevant. The other approach is real-time feedback 
control (RTFC) [391-396,638], in which the same quantum system is followed in real 
time around the feedback loop, and the measurement (or interaction with a quantum 
"controller" ) has a significant effect on the system's evolution. 

There are two distinct types of RTFC, which differ by the nature of the controller. 
In one approach to RTFC, measurements are employed to probe the quantum system, 
and the gathered information is processed classically off-line in real time to assess the 
best, next control action [391-396]. The evolution of the controlled quantum system is 
governed by two effects: coherent (unitary) action exerted by the classical controller^" 
and incoherent (non-unitary) back action exerted by the measurements. A generalized 
description of measurement-based RTFC employs quantum filtering theory [639-641]. 
Recently, another type of RTFC — referred to as coherent feedback control* — has 
drawn much attention [638, 642-648]. In this approach, no measurements with a 
classical output signal are performed; instead, an ancillary quantum system serves as 
the controller. The controller influences the evolution of the system of interest via a 
direct interaction between them. Additionally, external classical forces also can be used 
to act upon the system, the controller, or both. The system of interest together with 
the controller are characterized by the entirely quantum nature of the information flow 
— coherence is not destroyed by measurements, which is the source of the name given 
to this type of control [638,642]. Coherent feedback control can be viewed as a quantum 
analog of Watt's flyball governor — an automatic self-regulating quantum machine [649]. 

It was recently shown [650] that the evolution of a quantum system undergoing any 

+ The free evolution of the system can be included while the off-line modelling is performed. 

* We will also use the term coherent RTFC. The choice of terminology is standard in the field and 

should not be confused with the notion of coherent control employed in AFC and QOCT. 
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type of RTFC (or a combination thereof), including effects of measurements, feedback 
actions and interactions with auxihary quantum systems, can be generally represented 
by a Kraus map. This result has an important consequence for the optimization 
complexity of RTFC, since observable-control landscapes of open quantum systems with 
Kraus-map dynamics are free from local traps (under the controllability assumption) 
[273] . Remarkably, this trap-free landscape topology unifies virtually all types of control, 
including measurement-based RTFC, coherent RTFC, AFC, and open-loop control (the 
Kraus-map description of open-system dynamics is also generally valid for the latter 
two types when the controlled system is coupled to an environment). The possibility of 
employing this general unifying feature of controlled quantum dynamics for development 
of hybrid control schemes is discussed in section 8.4. 

Although the history of AFC in the quantum realm is short (it was proposed in 
1992 and first experimentally implemented in 1997), it has become a well established 
and popular laboratory tool successfully employed in growing numbers of experimental 
studies (see section 5). In contrast, while RTFC of quantum systems was first 
proposed in 1983 [391], its laboratory implementation [642,648,651-653] thus far has 
been much less extensive. Implementing RTFC on the atomic or molecular scale 
encounters significant technical difficulties. One especially important obstacle is that 
many interesting quantum phenomena occur on a time scale which is too short to allow 
for processing of the measurement data in classical controllers based on conventional 
electronics (i.e., the issue of loop latency). Coherent RTFC can overcome the latency 
issue, but in this case the controller itself may require precise engineering to assure 
quality control performance of the quantum system. Nevertheless, interest in potential 
applications of RTFC for manipulation of quantum systems is high and the theoretical 
activity in this field is growing [654-665]. Of particular importance (especially for 
applications in quantum information sciences) is the ability of RTFC to stabilize 
the dynamics of quantum systems in the presence of noise [646, 648, 666-669] and 
achieve robust control performance in the presence of uncertainties in the system 
Hamiltonian [670]. Other interesting possibilities include the use of RTFC for quantum 
error correction [671-676], generation and protection of entanglement [668,677-682], 
cooling of quantum systems [651,68-3-686], and quantum state purification [687-693]. 
In the related field of quantum metrology, real-time feedback was employed to approach 
fundamental quantum limits of measurement accuracy [694,695]. Both theoretical and 
experimental aspects of RTFC should continue to draw significant attention in the 
future. 

8. Future directions of quantum control 

Common sense dictates that the future is notoriously difficult to predict, but it is also 
the nature of science to try and anticipate new directions that will expand current 
knowledge. The evident paths followed in the development of the quantum control field 
during the last two decades provide a basis for projection, with due caution, on how 
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current trends may evolve upon going forward. Below we try to identify some critical 
theoretical and technological issues, where breakthroughs are required to significantly 
increase the capability of controlling quantum phenomena. 

8.1. Input-output maps for quantum control simulations 

Except for the special situation of measurement-based RTFC (where measurement back 
action is a distinctively non-classical feature), one may naively conclude that there 
are no fundamental differences between designing controls for quantum and classical 
systems. The distinctions seem to lie in solving classical equations of motion in one 
case and the Schrodinger equation in the other, but otherwise the method of finding the 
optimal control solution is basically the same. However, from a practical perspective, 
the difference between solving classical and quantum equations of motion is fundamental 
due to the exponential increase of the Hilbert space dimension characteristic of quantum 
systems. This is the reason why simulating controlled quantum dynamics of multi- 
particle systems is so difficult. 

A qualitative breakthrough in open-loop control of complex quantum systems would 
be possible, if a way could be found to replace the laborious calculations of quantum 
dynamics with "black-box" models that essentially capture the main features of the 
processes leading to a particular control objective (e.g., breaking of a specific molecular 
bond). The goal is to perform a modest number of simulations and use the information 
to generate an input-output map from the applied control field to the resultant change 
in the control objective. In this fashion, the input-output map aims to capture the 
relationship between the control and the system's reaction to it. This approach is 
commonly used in classical control problems in many areas of engineering; however, at 
the present time, we do not know how to effectively determine these input-output maps 
for complex quantum systems, such as molecules. 

An example of a method proposed for identifying nonlinear input-output maps 
for quantum control studies is the high- dimensional model representation (HDMR) 
technique [696-700]. The total number of points in the search space for a quantum 
control optimization problem (and for many other problems in science and technology) 
grows exponentially with the number of input variables (this situation is sometimes 
called the "curse of dimensionality"). In HDMR, the input-output map is characterized 
by a hierarchy of contributions from the input variables acting independently, in pairs, 
triples, etc. For many important problems, with an appropriate choice of the variables, 
only low order input variable cooperativity is significant. This property can be used 
to dramatically reduce the effort required to explore the map. Approaches such as 
HDMR are designed for systems with a large number of input variables with the aim of 
learning the input-output map using a number of simulations that grows relatively slowly 
(e.g., polynomially) with the number of input variables [700]. Specifically, the use of 
nonlinear functional HDMR-type maps for quantum control problems was discussed by 
Geremia et al. [701]. Such input-output maps would be of value as well when generated 
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from experimental data, as they would constitute the control landscape. Although 
much is now understood about the topology of quantum control landscapes, there is 
little information about non-critical point structural features. 

Another approach popular for constructing nonlinear input-output maps employs 
neural networks [702,703]. Recently, neural networks were used to model ultrafast laser 
control of SHG, molecular fluorescence yield, and photoelectron spectra from resonant 
strong-field ionization of potassium atoms [704,705]. However, it was found that the 
amount of data required for the training of a neural network significantly increases with 
the complexity of the correlations which are to be modelled. While reproduction of the 
training data worked very well, extrapolation to regions of the parameter space which 
were not covered by the training data was a challenge. 

Development of efficient and accurate input-output maps for control of complex 
quantum phenomena remains an important objective. Ideally, after a modest effort at 
learning a map, it could then be used in a highly efficient manner to seek out new 
controls and dynamical regimes with favorable characteristics. 

8.2. Analysis of quantum control landscapes 

The introduction of quantum control landscapes in the last few years is an important 
theoretical advance in the field. The nature of the control landscape topology has direct 
implications for the ease of finding effective controls in the laboratory. The analysis 
of the control landscape topology and other structural features can provide the basis 
for investigating the complexity of optimizing different types of control objectives. In 
turn, this understanding can help identify the most suitable optimization algorithms for 
various theoretical and experimental applications of quantum control (see sections 4.4, 
4.5, and 5.19). In addition, the landscape analysis may be extended to the study of 
quantum control problems involving simultaneous optimization of multiple objectives 
(see sections 4.3 and 4.4). This research area is still rapidly developing with much 
remaining for investigation. In particular, an open issue that deserves significant 
attention is the effect of field constraints (e.g., due to limited laboratory resources) 
upon the accessible regions of quantum control landscapes. 

There are several additional research directions for which the analysis of the control 
landscape features may provide important insights. One ubiquitous problem with wide- 
ranging implications is evaluation of the robustness of control solutions to noise and 
imperfections, which depends on the degree of fiatness of the control landscape around 
an optimal solution. Another interesting issue is related to a phenomenon discovered for 
observable control of an open quantum system prepared in a mixed state and coupled 
to a thermal environment [273]. Specifically, the range of the control landscape (i.e., 
the difference between the maximum and minimum expectation values of the target 
observable) decreases when the temperature of the environment raises. Therefore, 
an important application of control landscape analysis would be determination of the 
fundamental thermodynamic limits on the control yield for open quantum systems. 
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8.3. Future applications of AFC 

As discussed in section 5, AFC has proved to have broad practical success as a means 
for achieving optimal control of quantum phenomena in the laboratory. Particularly 
impressive is the breadth of applications, ranging from optical systems, to atoms, to 
semiconductor structures, to biologically relevant photochemical processes in complex 
molecules, etc. One clear trend is the extension of AFC applications towards the 
manipulation of increasingly more complex systems and phenomena. Along this avenue, 
implementation of AFC could bring significant benefits to such areas as near and even 
remote detection of chemical compounds (first steps in this direction have been recently 
made [490,706]), optimal control of molecular electronics devices, and optimal control 
of photochemical phenomena in live biological samples (including nonlinear microscopy 
and ODD, as indicated by several recent experiments [539-542]). 

We can also envision increasing use of AFC for optimal quantum control of 
photophysical phenomena. One important area is coherent manipulation of quantum 
processes in solid-state systems, especially in semiconductor quantum structures 
[707,708]. Another potential application is optimal storage and retrieval of photonic 
states in atomic- vapor and solid-state quantum memories [358,359,709-714]. The AFC 
methodology may be also applicable to physical problems where, instead of laser pulses, 
other means (e.g., voltages applied to an array of electrodes) are used to implement 
the control. Examples could include optimal control of coherent electron transport in 
semiconductors by means of adaptively shaped electrostatic potentials [715], coherent 
control of charge qubits in superconducting quantum devices by gate voltages [716,717], 
and coherent control of photonic qubits in integrated optical circuits via the thermo- 
optic effect [718]. Several types of quantum systems (e.g., flux qubits in superconducting 
quantum devices, hyperfine-level qubits in trapped neutral atoms and ions, electron 
spins of donor atoms in silicon, etc.) can be controlled by pulses of microwave radiation 
(e.g., AFC-optimized dynamical decoupling [594] of trapped- ion qubits by a sequence 
of microwave vr pulses was discussed in section 5.18). Many possible applications 
of AFC could have significant implications for the progress in the field of quantum 
information sciences. A new domain of quantum control involves manipulation of 
relativistic quantum dynamics with extremely intense laser fields [719] for accelerating 
particles and even intervening in nuclear processes in analogy with what is happening 
in atomic-scale control. It is reasonable to forecast that AFC methods could become 
useful for optimal control of such laser-driven high-energy phenomena. 

Despite significant advances achieved in the field of quantum control during the 
last decade, the experimental capabilities are limited by currently available laser 
resources. It is likely that existing practical limitations, in particular, the relatively 
narrow bandwidth of femtosecond lasers, restrict the achievable yields in some AFC 
experiments. One might expect that many new applications would open up, if reliable 
sources of coherent laser radiation with a much wider bandwidth became available. Such 
resources could make possible the simultaneous manipulation of rotational, vibrational. 
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and electronic processes in molecules in a more effective fashion, thereby achieving 
hitherto unattainable levels of control. Moreover, if pulse-shaping technology can be 
extended to coherent radiation in the attosecond regime as well as in the range of MeV 
photon energies, a multitude of new applications in X-ray spectroscopy, medical physics, 
and control of nuclear dynamics could arise. 

8.4- Hybrid methods of quantum feedback control 

Despite the significant technological difficulties on the path to routine practical 
application of RTFC, the potential benefits are alluring (see section 7). An interesting 
question is whether AFC, whose practical utility has already been well established, 
can be used to aid in the implementation of RTFC (measurement-based, coherent, or 
both). Due to the apparent technological differences between AFC and RTFC, thus 
far they have been considered as separate branches of quantum control. However, 
it has been recently shown [650] that AFC and RTFC share a common fundamental 
landscape topology characterized by the absence of local traps (i.e., all sub-optimal 
extrema are saddles provided that the controllability condition is satisfied). Since the 
control landscape topology strongly influences the optimization complexity, this finding 
may have immediate practical importance. Moreover, the unification of the seemingly 
different AFC and RTFC approaches at a fundamental level suggests the possibility of 
developing new laboratory realizations that combine these currently distinct techniques 
of quantum feedback control in a synergistic way. For example, some form of AFC 
might be used to optimize the design or construction of quantum controllers employed 
in coherent RTFC. Development of "hybrid" quantum control schemes incorporating 
both AFC and RTFC (in particular, for control and stabilization of quantum computing 
systems) could provide significantly enhanced flexibility in the laboratory. 

8.5. Material control 

In addition to the manipulation of quantum dynamics via application of optimal 
external fields, there is the prospect of performing material control through alteration 
of the internal system properties (i.e., the spatial structure or matrix elements of the 
system Hamiltonian) , with the aim of identifying optimal materials and system designs. 
Analogous to the circumstance of a particular quantum system acted upon by a family 
of homologous external control fields, we can consider the controlled response of a family 
of homologous quantum systems to a particular field. In the former case, a control level 
set consists of all homologous control fields that produce the same expectation value of 
the target observable when applied to a particular quantum system. This level set can 
be explored, for example, by homotopy trajectory methods (e.g., D-MORPH), in order 
to identify control solutions with desired properties (see section 4.4). In the latter case 
of material control, a level set consists of all dynamically homologous quantum systems 
that produce the same expectation value of the target observable when controlled by 
a particular field. For example, each quantum system may be specified by a point in 
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a hypercube whose edges are labeled by Hamiltonian matrix elements. A variation 
of the D-MORPH method can be used to explore a system level set by continuously 
warping the corresponding Hamiltonian [720]. At this juncture little is known about 
either homologous control fields or homologous quantum systems. Exploration of these 
topics could reveal the systematic aspects of control over quantum phenomena. 

Morphing through Hamiltonian structure in the laboratory can be physically 
realized in many different ways, with broad and yet largely unexplored possibilities. 
For example, the properties of light-sensitive materials could be varied using families 
of structurally similar chemical compounds, characteristics of doped semiconductors 
can be varied by changing the concentration of dopant atoms and the depth of 
implanting, etc. Material control is potentially applicable to a wide set of problems 
in various areas of science and technology. Possible applications include, for example, 
development of photodetectors with higher efficiency and faster response time, molecular 
switches with increased sensitivity and durability, quantum computing systems with 
enhanced immunity to environmentally-induced decoherence and improved robustness 
to instrumental noise, etc. Exploiting the accessible variations in Hamiltonian structure 
as a means for achieving optimal quantum control is a potentially important area of 
future research, including exploration of the associated control landscapes, development 
of adaptive and open-loop techniques, investigation of effective methods of combining 
material and electromagnetic control, and adaptation of the theoretical concepts to 
various practical applications. 

8. 6. Scientific and engineering goals of optimal quantum control 

The general goal of science is to understand nature, including the structure of physical 
systems and characteristics of the system dynamics, while the goal of engineering is 
to design and implement a system that will function in a prescribed manner in the 
best possible way. Optimal quantum control draws together science and engineering to 
incorporate both goals: (1) to understand the dynamical behavior of quantum systems 
and the mechanisms by which these processes can be managed, and (2) to require optimal 
functional performance through the achievement of prescribed control objectives in the 
best possible way. 

An important feature evident in the prior development of quantum control is the 
impact of progress in one aspect of the subject on advancing another. We expect 
that this trend will continue in the future, as a better understanding of the underlying 
physical processes would aid in choosing better control tools and thereby achieving 
a higher degree of performance. In turn, the ability to steer system evolution in an 
optimal fashion should facilitate the acquisition of knowledge about the underlying 
control mechanisms and other properties of the system. For example, in many AFC 
experiments, the characteristics of the resultant optimal control fields were used (often 
in combination with additional measurements and/or simulations) to decipher physical 
mechanisms responsible for the achieved control [459-463,465,467,468,472,474-480,501, 
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518,520,521,533,535,537,557,558,560,567-569]. Also, a method for analysis of quantum 
control mechanisms through Hamiltonian encoding was recently developed [721-728] and 
experimentally tested [720]. However, much additional theoretical and experimental 
work is still needed to better understand the controlled dynamics of complex systems, 
especially those in the condensed phases. 

Recent experiments [730-734] discovered manifestations of long-lived electronic 
quantum coherence in energy transfer processes in light-harvesting complexes of 
photosynthetic systems. Evidence of long-lived electronic and vibrational quantum 
coherence was also found in intrachain energy transfer in a conjugated polymer [735,736]. 
These findings raise an important question about the role of coherent quantum effects in 
the dynamics of energy transfer and other photoinduced processes in complex chemical 
and biological systems. This issue already attracted significant attention [737-740], 
and a number of theoretical models [741-761] were developed to explain the existence of 
quantum transport in the presence of strong coupling to a thermal environment. Further 
advances in this area may lead to a better understanding and more effective control of 
photophysical and photochemical quantum phenomena in the condensed phase. 

Hamiltonian identification is a potentially important application of quantum control 
aimed at revealing detailed information about physical systems. Extraction of the 
Hamiltonian from measured data is an inverse problem that generally suffers from 
being ill-posed (i.e., the Hamiltonian information is unstable against small changes of 
the data), which arises because the data used for inversion are inevitably incomplete. 
Recent attempts to address this challenging problem include Hamiltonian identification 
via inversion of time- dependent data (instead of traditional use of time-independent 
spectroscopic data) [762-767] and application of global, nonhnear, map-facilitated 
inversion procedures [768,769]. In this context, it appears that suitable controls can be 
used to significantly increase the information content of the measured data. For example, 
it may be possible to control the motion of a molecular wave packet to gain more 
information on interatomic forces in selected regions of a potential energy surface [14]. 
This concept has seen some nascent development by Geremia and Rabitz [770, 771] 
who proposed the notion of optimal Hamiltonian identification (OHI). OHI aims to 
employ coherent control of quantum dynamics to minimize the uncertainty in the 
extracted Hamiltonian despite data limitations such as finite resolution and noise. The 
proposed OHI implementation operates in a manner similar to an AFC experiment, 
using closed-loop optimization guided by a learning algorithm to discover controls that 
minimize the dispersion of the distribution of Hamiltonians consistent with the measured 
data. Numerical simulations indicate that an optimal experiment can act as a tailored 
filter to prevent laboratory noise from significantly propagating into the extracted 
Hamiltonian [771]. Ideally, upon each cycle of the experiment more information about 
the Hamiltonian will be extracted, which in turn would be used to guide the next cycle, 
etc. A critical component of OHI is the need for real-time numerical simulations of the 
quantum system's dynamics on the fly. OHI will require further development [203] of 
inversion algorithms, computational capabilities (e.g., possibly including input-output 
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maps discussed in section 8.1), and experimental techniques to become practical. 
9. Concluding remarks 

It would be impossible to cover in a paper of any reasonable length all of the advances 
that have been made in the last two decades in the field of quantum control. Some areas 
that did not receive full attention here were considered in more detail in other review 
articles and books (in particular, those cited in section 1), to which we refer the interested 
reader. For example, thematic reviews are available on control via two-pathway quantum 
interference [1,8,21,24], pump-dump control [7,11], control via STIRAP [10,16], control 
via WPI [42], the formalism and applications of QOCT [14,36,38,40], controllability 
of quantum systems [36], the formahsm of quantum control landscape theory [39], 
femtosecond pulse-shaping technology [4,5, 12,27], femtosecond laser control of X-ray 
generation [35,41], quantum control experiments with "rational" control designs [30], 
quantum control applications in nonlinear spectroscopy and microscopy [34,44], and 
control of quantum dynamics on the attosecond time scale [45]. While we tried to 
provide a comprehensive account of laboratory AFC of quantum phenomena, more 
detailed discussions of some important AFC experiments are available in earlier reviews 
[15,18,19,25,26,31,34,35,37,41]. New papers, often containing significant theoretical 
and experimental results in quantum control, appear now almost on a daily basis. 

In this paper, our goal was to give a perspective and prospective on the field 
highlighting the evolution of important trends in quantum control. A look into the past 
together with a review of current, cutting-edge research were used to cautiously forecast 
topics of future interest. We also attempted to emphasize the synergistic connection 
between the theoretical and experimental advances, which has been immensely beneficial 
for the development of the field. We believe that sustaining this productive interplay 
between theory and experiment will be critical for future progress. This paper aimed to 
provide the basis to better understand which aspects of theoretical research are having 
a high impact on laboratory control of quantum phenomena. At the same time, a 
complementary goal of this work was to point out the experimental aspects of quantum 
control that have special significance and relation to theory. Although the scope of 
experimental and theoretical research in quantum control is vast, we hope that this 
work provides a valuable bridge for the community involved as well as for those outside 
who are interested in understanding the reasons for the fervor in the field. 
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